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ABSTRACT 


"  The  purpose  of  this  project  was  to  determine  the  temperatures 
and  luminosities  of  the  individual  components  of  eclipsing  binary 

star  systems.  ~^r.  Richard  L.  Walker  of  the  U.S.  Naval  Observatory 

J 

provided  a  list'  of  such  systems  which  were  as  yet  undetermined. 

The  information  sought  was  gained  by  UBV  photometry  of  a  system  at 
total  eclipse  and  at  a  time  outside  eclipse.  The  light  at 
totality  is  due  entirely  to  the  occulting  star,  and  outside 
eclipse,  both  stars  contribute  fully.  A  method  is  derived  for 
subtracting  out  the  light  of  the  occulting  star  to  obtain  measure¬ 
ments  of  the  occulted. 

A  recently  published  technique  by  William  E.  Harris  was  used 

/ 

to  reduce  raw  observational  data.  Essentially  a  simultaneous, 

/ 

multi-linear  solution  for  both  atmospheric  extinction  and  instru¬ 
mental  transformation,  this  method  was  programmed  onto  the  Naval 
Academy's  extensive  computing  facilities.  It  was  found  that  this 
method,  as  claimed,  gives  accurate  results  with  a  minimum  of  time 
and  effort,  upon  completion  of  the  initial  phase  of  programming. 

^Systems  for  which  a  complete  solution  (temperature  and  lumi¬ 
nosity  of  both  components)  was  reached  include:  TU  Camelopard!, 

TW  Draconis,  AK  Herculis,  V566  Ophiuchi,  W  Ursae  Majorls,  and  AG 
Virginia.  Systems  observed  only  during  totality,  thus  solving 
only  the  occulting  star,  include  a  Corona  Borealis  and  AM  Leonis. 
RS  Canes  Venatici  and  TZ  Bootes  were  observed  only  out  of  eclipse, 
and  must  await  further  study.  Once  a  solution  for  a  system  was 
obtained,  it  was  presented  graphically  on  a  Hertzsprung-Russell 
diagram,  and  was  examined  from  the  viewpoint  of  binary  evolution. 
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I  -  THEORETICAL  CONSIDERATIONS 

A  primary  goal  of  binary  scar  astronomy  Is  Co  understand  the 
circumstances  and  mechanisms  which  led  to  their  formation  and 
which  will  dictate  their  future.  Definitive  solutions  to  these 
problems  clearly  do  not  yet  exist,  and  are  for  the  most  part 
beyond  the  scope  of  this  paper.  Yet,  to  make  any  real  progress  in 
this  direction,  one  requires  a  firm  grasp  of  the  nature  of  binary 
stars.  This  section  will  therefore  discuss  the  dynamics,  physical 
characteristics',  and  evolution  of  binary  systems,  prefaced  by  a 
treatment  of  similar  topics  for  ordinary  single  stars,  since  these 
form  the  components  of  binary  systems.  These  theoretical  con¬ 
siderations  will  then  lay  the  foundation  for  interpretation  of 

r 

experimental  results  and  hopefully  lead  to  insights  into  binary 
system  development. 


Single  Stars 

For  the  purposes  of  this  paper,  stars  possess  two  essential 
properties.  The  first  of  these  is  temperature.  Unfortunately, 
the  most  crucial  temperature,  that  which  controls  the  type  and 
rate  of  energy  generation,  remains  unobservable  in  the  central 
regions  of  the  star.  On  Earth,  one  can  measure  only  the  tem¬ 
perature  at  the  depth  in  from  the  surface  at  which  the  star  becomes 
opaque  to  the  wavelength  of  light  observed.  If  we  limit  the 
range  of  wavelengths  to  the  visible  and  the  near  ultra-violet 
regions,  the  differences  in  penetration  depths  relative  to  stellar 
radii  becomes  insignificant,  and  we  have  defined  a  useful  measure 
of  temperature  with  which  to  compare  stars,  the  effective  temperature. 
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There  are  several'  alternative  measures  of  temperature. 
Spectral  class  orders  scars  by  the  characteristics  of  their 
spectra,  which  depend  in  turn  on  the  excitation  levels  of  the 
various  atoms  in  their  atmosphere.  This  is  a  standard  way  to  spe¬ 
cify  temperature,  primarily  for  economy  of  print.  The  spectral 
classes  are  described  in  table  1-1. 


TABLE  1-1 


SPECTRAL  CLASS 

T(x103K) 

CHARACTERISTICS 

05 

50.0 

Absorption:  iodized  He  ,  N,  0,  Si 

BO 

27.0 

Absprption:  He 

B5 

16.0 

He4  declining,  Balmer  series  Increasing 

AO 

10.4 

No  He  lines,  H  lines  dominant 

A5 

8.2 

Singly  ionized  Ca  increasing 

F0 

7.2 

U  weakening,  Ca  increasing 

F5 

6.7 

Neutral  metal  lines  increasing 

GO 

6.0 

Ca  strong,  many  neutral  metals 

G5 

5.5 

CN,  CH  molecular  bands 

K0 

5.1 

Ca  maximum,  neutral  metals  stronger 

K5 

4.3 

Many  molecules 

MO 

3.7 

Neutral  metals  intense,  U  weak 

M5 

3.0 

TiO  appears 

Color  also  measures  stellar  temperatures,  ranging  from 
bluish-white  for  hot  0  and  B  stars  to  yellow  for  G  stars  (like  the 
Sun)  to  red  for  cool  M  stars.  Unfortunately,  since  determining 
the  color  of  an  unknown  star  is  largely  subjective,  it  is  not  a 
strictly  precise  temperature  measurement.  However,  color  index, 
the  difference  in  intensities  of  starlight  at  two  specified  wave¬ 
lengths,  does  give  an  accurate  result  to  the  approximation  that 
stars  radiate  as  black  bodies.  One  color  index  uniquely  deter¬ 
mines  a  Planck  curve,  which  specifies  the  temperature.  Addition 
of  a  third  wavelength  (a  second  color  index)  shows  that  a  star  is 


not  a  perfect  black  body,  and  reveals  subtleties  of  the  stellar 
atmosphere,  which  will  later  become  quite  Important. 

The  second  essential  property  of  a  star  is  its  luminosity, 
the  total  energy  emitted  per  unit  time: 

L-(4*r2)(oT4)  (1-1) 

where  o  is  the  Stef an-Boltzmann  constant.  Since  stellar  radii  are 
generally  not  known,  this  equation  is  of  limited  usefulness  in 
finding  luminosities. 

The  concept  of  stellar  magnitudes  allows  one  to  obtain  a 
usable  measure  of  stellar  luminosities.  The  magnitude  scale  ranks 
stars  on  the  basis  of  visual  response  to  their  relative  brightnesses. 
Due  to  ocular  physiology,  this  scale  is  logarithmic: 

m2~mi*2.51og(Li/L2)  (1-2) 

We  notice  from  this  equation  the  unfortunate  historic  convention 
that  brighter  stars  have  lower  magnitudes.  A  first  magnitude  star 
is  one  hundred  times  brighter  than  a  sixth  magnitude  star. 

Since  the  flux  of  stellar  radiation  obeys  an  Inverse  square 
law,  stars  closer  to  the  Earth  will  appear  unnaturally  brighter 
than  those  further  away.  To  equalize  this  distance  dependency,  we 
define  absolute  magnitude  as  the  magnitude  of  a  scar  if  it  were 
ten  parsecs  (about  thirty-two  light  years)  distant.  As  an 
example,  the  proximity  of  Che  Sun  gives  it  an  apparent  magnitude 
of  -26.8,  whereas  its  absolute  magnitude  is  only  +4.8. 

Once  we  know  a  star's  temperature  and  luminosity,  we  can 
locate  a  star's  position  on  a  Hertzsprung-Russell  diagram,  a  plot 
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of  luminosity  versus  temperature,  as  in  figure  l-l1 .  An  Important 
feature  of  this  plot  is  that  stars  are  not  randomly  distributed  as 
one  might  expect,  but  rather  group  themselves  along  specific 
bands.  These  bands  are  termed  luminosity  classes,  and  are  the 
subtle  distinctions  measured  by  a  second  color  index. 

The  B-R  diagram  also  gives  a  graphical  representaiton  of  a 
star's  evolution.  At  any  time  in  its  life,  a  star  will  occupy  a 
specific  position  on  the  diagram.  If  we  connect  all  the  points  of 
its  lifetime,  we  obtain  an  evolutionary  track.  The  next  task  is 
thus  to  approach  these  tracks  theoretically,  searching  for  the 
signiflgance  of  and  explanations  for  the  grouping  found  on  the  B-R 
diagram. 


Pigure  1-1. 
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Ic  is  generally  accepted  chat  stars  originate  in  areas  where 
the  local  density  of  interstellar  gas  and  dust  is  high  enough  to 
initiate  gravitational  condensation.  As  this  cloud  contracts,  its 
particles  lose  gravitational  potential  energy.  The  virial  theorem 
for  inverse  square  interactions: 

<T>— V2<U>  (1-3) 

shows  that  half  of  this  energy  loss  goes  into  increasing  the 
average  kinetic  energy  (i.e.  the  temperature)  of  the  cloud.  The 
other  half  must  be  radiated  away.  Initially,  the  mean  free  path 
of  these  photons  is  large,  and  they  can  escape  from  the  cloud 

regardless  of  their  depth  of  emission.  This  contraction  stage  is 

••  3  2 

relatively  rapid,  on  the  order  of  10  years  ,  less  for  more 

massive  clouds  and  greater  for  lighter  clouds.. 

At  this  time,  however,  Che  density  of  the  proto-star  has 
increased,  lowering  the  mean  free  path  of  the  photons  such  that 
they  will  most  probably  Interact  and  release  their  energy  well 
inside  the  proco-star.  This  radiation  pressure  increases  the 
internal  temperature  and  establishes  a  large  inwardly  increasing 
density  gradient.  The  contraction  also  slows  due  to  Che  radiation 
pressure  until  the  internal  temperature  reaches  1300K,  at  which 
point  H2  dissociates3  (endothermically),  again  allowing  a  rapid 
contraction.  This  phase  continues  until  a  state  of  approximate 
hydrostatic  equilibrium  is  reached. 

After  about  107  years,  nuclear  reactions  begin  to  contribute 
to  the  proto-star's  luminosity.4  These  reactions  are  only  precur¬ 
sors  to  the  eventu <1  energy  reduction  mechanisms,  since  internal 
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temperatures  are  still  too  low  to  ignite  complete  fusion  chains. 
Figure  1-2S  gives  a  composite  view  of  stellar  life  for  masses  from 
one-half  to  fifteen  solar  massses  to  this  point. 

Before  following  stellar  evolution  further,  it  is  best  to 
momentarily  digress  to  the  subject  of  degenerate  matter  and  its 
effects  on  the  subsequent  life  of  stars.  For  an  n-dimenaional 
object,  there  exists  a  2n-dimen.aion.al  phase  space  whose  co¬ 
ordinates  are  the  a  spatial  co-ordinates  (q^)  and  the  respective 
momenta  co-ordinates  (p^).  The  Heisenberg  uncertainty  rela¬ 
tionship  between  position  and  momentum  gives  the  mimlmun  volume  of 
a  state  in  2n  space: 


(nAq1)(IIhp1)>(h/  2ir)n  (1-4) 

into  which  no  two  fermions  (spin  one-half  particles,  for  our  pur¬ 
poses,  electrons  and  neutrons)  can  be  placed,  due  to  Pauli  exclu¬ 
sion.  The  available  physical  (spatial)  volume  is  controlled  by 
the  particle  density;  an  Increase  in  density  causes  the  average 
particle  separation  to  decrease  (decreasing  HAq^),  limiting  the 
number  of  available  states.  Likewise,  volume  in  momentum  space  is 
controlled  by  temperature;  an  increase  in  temperature  raises  the 
average  particular  momentum,  thereby  Increasing  HAp^.  When  all 
the  available  states  for  a  given  density  and  temperature  are 
occupied,  we  have  degenerate  matter,  and  the  particles  exert  a 
great  pressure  (degeneracy  pressure)  to  avoid  violating  the  Pauli 
exclusion  principle. 

Thus,  a  star  isothermally  contracting  will  eventually  stop; 
the  number  of  momentum  states  are  constant,  while  the  number  of 
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spatial  states  decreases.  The  total  number  of  states  therefore 
decreases,  while  the  number  of  particles  remains  the  same,  making 
electron  degeneracy  inevitable.  However,  if  we  increase  the  tem¬ 
perature,  more  momentum  states  become  available,  allowing  the 
number  of  spatial  states  to  continue  to  decrease.  Thus  we  get  the 
curious  result  that  adding  energy  to  degenerate  matter  causes  a 
contraction,  as  long  as  degeneracy  is  not  removed. 

The  proto-star  we  left  in  figure  1-2  now  has  two  mass-dependent 
alternatives.  If  its  mass  is  less  than  a  tenth  of  a  solar  mass, 
the  internal  temperature  never  reaches  the  ignition  point  of 
hydrogen  fusion.  This  "black  dwarf"  releases  energy  through  the 
motion  of  its  particles,  gradually  cooling  to  equilibrium  with  the 
space  around  it.  If  its  mass  is  greater  than  this  limit,  the  star 
will  commence  its  long  life  of  hydrogen  burning. 

The  onset  of  hydrogen  fusion  is  termed  Zero  Age  Main 

Sequence.  (The  main  sequence  is  the  densely  populated  diagonal 

band  on  the  H-R  diagram.  Most  of  the  stars  we  see  are  on  the  main 

sequence  simply  because  stars  spend  most  of  their  lives  burning 

hydrogen.)  At  this  point,  the  chemical  composition  of  the  star 

can  be  considered  homogeneous,  with  typical  values  of  7 1%  H,  27% 

He,  and  2%  heavier  elements,  by  weight.6  The  star  generates 

1  *♦ 

energy  by  fusing  four  H  into  one  He  ,  primarily  by  either  the 
proton-proton  chain  or  the  carbon-nitrogen  cycle  (see  Appendix  I). 
The  important  feature  of  both  reactions  is  that  only  0.71%  of 
the  mass  of  the  four  hydrogens  is  converted  into  energy.  Thus, 
during  the  main  sequence,  the  stellar  mass  remains  essentially 
constant,  but  the  composition  changes  drastically. 
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1  4 

As  4U  ♦He  ,  che  temperature  decreases  slightly  due  to  che 
increased  average  molecular  weight  (u)  of  the  gas  as  a  whole,  and  the 

luminosity  increases  for  the  same  reason.  (L  ranges  in  propor- 

7  5  4  7 

tionality  from  u  *  in  small  scars  to  u  in  large  stars.  )  The 

star  thus  moves  upward  and  to  the  right  a  small  amount  while  on 
the  main  sequence.  The  stellar  mass  determines  che  length  of  the 
star's  stay  here;  more  massive  stars  have  higher  internal  tem¬ 
peratures,  which  causes  a  faster  nuclear  reaction  rate  and  a  more 

rapid  hydrogen  depletion.  For  the  sun,  this  time  would  be  about 

9  7  8 

9x10  years,  and  for  a  star  five  times  as  massive,  3x10  years. 

As  time  passes,  hydrogen  becomes  exhausted  in  the  central 
regions,  leaving  an  inert  helium  core  and  a  hydrogen  burning  shell 
around  it.  As  more  hydrogen  is  depleted,  the  helium  core  grows 
and  eventually  begins  to  gravitationally  contract.  This  releases 
gravitational  potential  energy,  which  raises  the  temperature  and 
che  fusion  rate  of  the  burning  shell.  The  increase  in  energy  flux 
from  this  shell  pushes  the  outer  regions  of  che  star  away  and  forces 
the  core  in  on  itself.  The  surface  temperature  decreases  with 
expansion,  but  the  luminosity  remains  more  or  less  constant  due  to 
the  increase  in  the  radiating  area. 

Meanwhile,  the  temperacure  and  density  of  the  contracting 
helium  core  continue  to  increase.  It  now  has  three  mass  and  tem¬ 
perature  dependent  alternatives.  If  the  star  has  less  chan  half  a 
solar  mass,  che  core  will  be  coo  cool  to  ignite  helium  burning, 

9 

and  will  suffer  electron  degeneracy.  The  outer  layers  will  con¬ 
tinue  to  expand,  leaving  behind  a  small,  degenerate  white  dwarf. 

This  releases  energy  through  che  motion  of  che  non-degenerate 


nuclei  (mn=> iOOOme ) ,  and  also  cools  co  equilibrium  with  its  surroundings. 

If  Che  scar  has  more  chan  chree  solar  masses,  cencral  tem- 
peracures  will  be  high  enough  Co  commence  helium  burning  via  the 
triple-alpha  chain  (see  Appendix  I).  This  new  source  of  energy 
expands  Che  core,  which  decreases  ics  oucward  pressure,  allowing 
Che  oucer  layers  Co  re-concracc.  A  rule  of  chumb  for  decermining 
qualitatively  Che  radius  variation  is  thac  each  nuclear  burning 
shell  reverses  Che  direction  of  expansion  or  contraction,  as 
illustrated  in  Figure  1-3. 10  The  temperature  in  Che  core  peaks 
and  then  decreases  with  expansion,  and  che  luminosity  decreases 
due  Co  Che  blanketing  effect  of  che  cooler  oucer  envelope.  Yec, 
as  che  envelope  contracts,  it  also  heats,  causing  che  luminosity 
Co  swing  back  up.  Just  as  in  che  hydrogen  burning  phase,  che 
helium  eventually  becomes  depleted,  leaving  a  predominantly 
carbon-oxygen  inert  core  surrounded  by  a  helium  burning  shell. 

This  iLOves  che  peak  temperacure  radially  outward,  further 
Increasing  che  luminosity. 

Scars  between  chese  two  mass  limits  ( V2m0<m*<  reach  che 
state  previously  described  for  massive  scars,  but  through  an 
entirely  different  mechanism.  At  che  termination  of  hydrogen 
burning  in  these  scars,  che  temperacure  is  sufficiently  low  to 
allow  Che  very  cencral  regions  Co  become  degenerate.  The  hydrogen 
burning  shell  pushes  the  core  in  on  itself,  starting  minor  helium 
burning.  This  added  energy  contracts  Che  degenerate  region, 
increasing  its  Cemperature  and  Chus  increasing  che  helium  fusion  rate. 
This  cycle  clearly  establishes  a  volatile  instability  as  long  as 
some  of  the  core  remains  degenerate.  This  is  termed  helium  flash; 
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Che  power  per  unit  mass  released  by  such  a  core  can  reach  10 
erg/gm/s.  The  sun  in  contrast  puts  out  10  erg/gm/s.11  The  tem¬ 
perature  soon  increases  enough  to  remove  degeneracy,  and  the  core 
expands  and  cools.  The  helium  core  now  burns  as  does  a  more 
massive  star's,  leaving  a  carbon-oxygen  inert  core  and  a  helium 
burning  shell. 

The  star  is  now  faced  with  three  further  mass  dependent 

alternatives:  a  degenerate  carbon-oxygen  white  dwarf,  ignition  of 

carbon  burning,  or  carbon  flash.  In  these  stages,  the  energy  loss 

through  neutrinos  has  an  appreciable  affect  on  the  critical 

masses;  at  least  six  solar  masses  are  required  to  avoid  carbon 

flash  ignoring  neutrino  loss,  but  fifteen  are  required  when  it  is 
12 

included.  Figure  1-4  summarizes  the  evolutionary  stages  related 
to  this  study.  Figure  1-513  shows  the  evolutionary  tracks  of 
stars  of  various  masses  from  Zero  Age  Main  Sequence. 

We  can  see  from  Figure  1-5  that  stars  of  moderate  masses  all 
move  into  the  giant  region,  and  very  massive  stars  become 
supergiants.  For  stars  less  than  three  solar  masses,  we  see  the 
violent  reaction  to  the  helium  flash  from  points  five  to  six. 

Each  doubling  back  of  a  star's  track  corresponds  to  the  ignition 
of  a  heavier  nuclear  fuel.  More  importantly,  we  can  see  that  the 
various  luminosity  classes  are  nothing  more  than  various  stages  of 
a  star's  life,  and  that  the  stars  above  the  main  sequence  all  have 
distended  atmospheres,  l.e.  high  volumes  at  low  densities.  This 
feature  becomes  critical  when  another  star  is  in  its  proximity,  as 
is  the  case  in  a  close  binary  pair. 
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Binary  Systems 

We  may  define  a  binary  system  as  two  stars  bound  In  orbital 
stability  by  their  mutual  gravitation.  Binaries  which  are 
resolvable  as  two  distinct  stars  In  a  telescope  are  termed  visual 
binaries.  Systems  which  reveal  themselves  by  Doppler  shifted 
spectra  (due  to  the  radial  component  of  their  orbital  velocity) 
are  spectroscopic  binaries.  Eclipsing  binaries,  the  focus  of  this 
paper,  betray  their  presence  by  periodic  changes  In  their  observed 
light  intensity.  The  eclipsing  system  holds  a  special  place  of 
Importance  In  astronomy,  since  they  are  the  most  easily  discov¬ 
erable  binary  type,  and  binaries  are  the  principal  source  of 
stellar  masses  .and  absolute  dimensions. 

Continuing  with  basic  definitions,  total  eclipse  occurs  when 
the  larger  star  completely  occults  the  smaller.  The  opposite  case 
is  annular  eclipse.  The  primary  eclipse  is  the  eclipse  of 
greatest  reduction  in  magnitude,  and  the  other,  shallower,  eclipse 
is  the  secondary.  The  light  curve  is  a  plot  of  intensity  against 
time,  which  shows  the  periodic  nature  of  the  alternating  eclipses, 
and  clearly  separates  primary  from  secondary  eclipse  based  on  the 
depth  of  the  minimum. 

Three  classes  of  eclipsing  binaries  can  be  defined  based  on 
characteristics  of  the  light  curve.  Detached  pairs  have  fairly 
constant  light  outside  the  eclipses  and  minima  of  usually  unequal 
depth.  Their  periods  are  on  the  order  of  several  days  to  years. 
Semi-detached  systems  have  smoother  light  curves  due  to  interac¬ 
tions  out  of  eclipse,  and  also  have  shallower  minima.  Their 
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periods  are  on  Che  order  of  days.  Contact  binaries  have  very 

shallow  ny<  rtimj)  of  almost  equal  depth,  and  eclipses  which  cover  a 

sizable  fraction  of  the  period,  which  is  on  the  order  of  a  few 
14 

hours.  Figure  1-6  shows  the  relative  dimensions  of  these  three 
types  of  binaries. 

Our  task  now  turns  to  understanding  the  mechanisms  which  vary 
the  light  outside  eclipses  in  semi-detached  and  contact  systems. 
Physical  distortion  is  the  simplest  of  these,  and  occurs  both  as 
tidal  distortion  due  to  mutual  gravitation,  and  as  equatorial 
distortion  due  to  the  rotation  of  the  star  Itself.  The  strong 
tidal  interaction  causes  both  deformation  and  synchronous  orbit,  so 
that  the  cross-section  area  radiating  towards  Earth  oscillates 

r 

with  the  orbital  period.  Since  the  tidal  force  elongates  the 
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Figure  1-6. 
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stars  along  Che  line  joining  their  centers,  chis  area  will  be  a 
minimum  at  eclipse  and  a  maximum  a  quarter  phase  lacer.  This 
tends  Co  distort  the  entire  light  curve  vertically.  The  rota¬ 
tional  distortion  has  an  effect  only  when  the  equatorial  planes  of 
the  Individual  components  are  Inclined  relative  to  the  orbital 
plane,  in  which  case  the  magnitude  of  this  effect  should  oscillate 
in  sign  with  the  system's  precessional  period15,  and  is  negligible 
on  our  time  scale.  Parallel  rotational  and  orbital  axes  cause  the 
distorted  rotating  star  to  always  present  the  same  cross-section. 

The  other  major  effect  on  light  outside  eclipse  is  the  irrad¬ 
iation  of  one  star  by  the  other,  sometimes  inaccurately  called 
Che  reflection  effect.  Radiation  falling  on  one  of  the  stars  can 
be  either  scattered  or  re-emitted.  In  the  first  case,  the  tem¬ 
perature  of  the  scatterer  is  not  appreciably  raised,  and  the 
resultant  radiation  is  partially  polarized.  In  the  second  case, 
the  temperature  of  the  stellar  atmosphere  is  raised  as  the 
Incoming  photons  excite  Che  atoms  to  various  higher  energy  levels. 
When  these  atoms  fall  back  down  to  lower  states,  they  emit  new 
photons,  but  not  necessarily  with  the  same  spectral  distribution 
as  before.  Both  these  effects  are  far  too  minor  to  affect  the 
internal  state  of  the  star. 

The  preceedlng  effects  alter  the  light  curve  primarily  out 
of  eclipse;  let  us  now  examine  variations  in  Che  eclipses  them¬ 
selves.  Limb  darkening  is  the  most  prominent  of  these.  This 
causes  the  parts  of  Che  surface  further  away  from  the  observer  to 
appear  darker  since  the  stellar  atmosphere  grows  opaque  at  a 
higher  (cooler)  altitude.  This  also  has  a  color  dependency,  as 
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optical  depth  is  a  function  of  wavelength.  Limb  darkening 
increases  the  initial  rate  of  intensity  decrease  during  eclipse, 
and  can  suggest  a  rule  for  determining  which  eclipse  is  total. 
During  totality,  we  see  only  the  occulting  star  (constant  light) 
until  the  smaller  star  re-appears;  but  during  an  annular  eclipse, 
we  see  the  smaller  move  across  the  larger.  Since  limb  darkening 
makes  the  center  of  the  larger  star  appear  brighter  than  its  edge, 
we  should  find  a  slight  ”11"  during  the  minimum  of  the  annular 
eclipse.  Thus  totality  occurs  when  the  minimum  is  flat. 

Gravity  darkening  affects  distorted  stars.  This  varies  the 
emergent  intensity  proportionally  to  the  local  gravitational 
acceleration16.  This  acts  independently  of  limb  darkening;  those 
regions  closer  to  the  volumetric  center  of  the  star  will  appear 
darker  than  those  further  away.  The  orientation  of  the  star  can 
thus  make  these  two  darkening  mechanisms  either  reinforce  or  can¬ 
cel  each  other's  contribution. 

To  obtain  a  more  useful  classification  of  binary  systems,  one 
which  is  also  better  suited  for  the  discussion  of  binary  evolution, 
we  must  make  use  of  the  Roche  model17.  We  envision  a  co-ordinate 
frame  whose  origin  is  at  the  system's  center  of  mass  and  which 
rotates  such  that  the  X-axis  always  joins  the  centers  of  the  two 
stars.  The  potential  U  at  point  (x,y)  becomes: 


U 


(1-5) 


where:  M  -  the  total  mass  of  the  two  stars 

f  “  the  fractonal  mass,  m^/m2  such  that  f< 1 
a  •  the  separation  of  the  two  stars 
ri«  the  distance  from  atari 
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Setting  U  equal  to  a  constant  and  solving  for  x  and  y  generates 

equipotential  surfaces  as  illustrated  in  Figure  1-7.  For 

large  U,  the  equipotentials  are  closed  ovals  about  the  masses.  As 

U  decreases,  the  ovals  grow  larger  until  some  Ucric  where  they 

touch  at  one  point,  L^.  We  call  this  surface  the  Roche  limit  or 

lobe.  Matter  at  point  will  be  in  unstable  equilibrium.  For  U 

less  than  Ucr^t,  both  stars  are  contained  within  the  same 

equipotential,  and  they  lose  their  binary  nature.  We  may  then 

think  of  the  Roche  limit  as  the  largest  closed  equipotential 

capable  of  containing  exactly  the  whole  of  each  component,  as 

1  8 

governed  by  the  mass  ratio,  f. 

We  can  now  use  the  Roche  limit  as  a  basis  for  defining  binary 
types.  A  detached  pair  has  both  components  well  within  their 
respective  lobes.  A  semi-detached  pair  has  one  star  filling  Its 
lobe,  while  the  other  is  within  its  limit.  A  contact  binary  has 
both  stars  filling  their  lobes,  and  can  be  considered  in  contact 
at  L}.  The  shape  of  the  Roche  lobes  suggests  great  tidal  distor¬ 
tion,  but  we  must  realize  that  only  the  outer  atmosphere  is 
greatly  distorted;  the  compact  nuclear  burning  central  region 
remains  unaffected,  and  evolves  as  only  a  function  of  the  star's 
mass. 

To  see  the  importance  of  Roche  lobes  to  the  evolution  of  a 
close  binary  system,  consider  what  happens  to  an  aging  star  whose 
outer  layers  are  distending.  Since  the  size  of  the  lobe  is  deter¬ 
mined  by  the  mass  ratio  and  the  separation,  both  of  which  have  no 
reason  to  change  while  the  stars  are  within  their  lobes,  the 
volume  available  for  expansion  is  constant  and  limited.  Thus,  the 
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outer  layers  of  one  star  will  eventually  distend  through  L} ,  and 

mass  will  transfer  between  the  two  stars.  This  will  in  turn 

change  the  mass  ratio,  and  thus  the  size  of  the  lobes.  This  effect 

1 9 

is  illustrated  (with  exaggeration)  in  Figure  1-8. 

Kopal  presents  a  detailed  account  of  the  mechanisms  of  mass 
20 

transfer.  Bypassing  the  mathematic  formalism,  he  predicts  dif¬ 
ferent  trajectories  for  a  particle  escaping  from  a  star  at  its 
Roche  limit  depending  on  the  particle's  ejection  velocity.  The 
particle  can  go  to  the  other  star,  directly  or  after  orbiting  the 
system;  it  can  return  to  the  original  star  after  passing  around 
the  other  star  one  or  more  times;  or  it  can  escape  the  system  alto¬ 
gether.  The  greater  the  ejection  velocity,  the  more  orbits  the 
particle  will  make. 

Some  of  these  paths  end  on  the  star  already  at  its  Roche 
limit.  This  star  cannot  stably  accept  additional  matter,  thus  the 
flow  will  increase  so  that  the  star  remains  at  its  limit.  Along 
with  an  increase  in  flow  rate  comes  increased  ejection  velocities, 
which  leads  to  more  particles  in  motion  around  or  escaping  from 
the  system.  This  model  therefore  predicts  clrcumstellar  activity 
in  the  form  of  streams  or  rotating  disks.  This  adds  another 
complication  to  the  light  received  on  Earth. 

Turning  to  the  question  of  binary  evolution,  we  find  three 
basic  theories  which  attempt  to  explain  binary  formation.  The 
first  proposes  that  two  initially  separate  stars  gravitationally 
capture  each  other.  However,  regions  of  high  stellar  density, 
where  such  a  capture  should  be  more  probable,  show  no  signs  of  an 
inordinately  high  percentage  of  binary  systems.  The  second  calls 
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for  che  fission  of  a  stellar  body  into  multiple  components,  but 

strong  density  gradients  in  internal  regions  of  star  prohibits 
21 

such  fragmentation. 

The  third  theorizes  a  pre-stellar  fission;  the  condensation 
of  a  contracting  cloud  about  multiple  mass  centers.  For  this  gra¬ 
vitational  instability  to  occur,  the  pre-stellar  cloud  must  have 
more  chan  one-hundred  solar  masses,  and  a  volume  of  about  a  cubic 
parsec.  We  can  see  why  there  should  be  fragmentation  by  examining 
the  cloud's  angular  momentum.  Due  to  both  the  large  radius  and 
the  effects  of  differential  rotation,  the  collapsing  cloud  has 
very  high  angular  momentum,  (on  the  order  of  107  times  that  of  a 
normal  star),  which  limits  its  contraction  as  a  single  body.  Yet 
the  cloud  wants  to  reduce  its  gravitational  potential  energy  as 
much  as  possible;  to  do  this  while  conserving  angular  momentum,  it 
may  fragment  into  many  sub -systems.  The  angular  momentum  is  con¬ 
served  by  the  revolution  of  the  fragments  about  the  total  center 
of  mass,  and  each  fragment  can  individually  contract.  We  are  left 

with  a  large  system,  mostly  in  two  to  four  star  sub-systems  with 
22 

stray  singlets.  All  these  stars  also  share  similar  initial  com¬ 
position  since  they  have  a  common  spawning  ground. 

Let  us  now  follow  the  evolution  of  a  hypothetical  binary, 
starting  after  both  components  have  started  nuclear  burning,  and 
are  both  within  their  Roche  lobes.  The  more  massive  star  will 
deplete  its  hydrogen  and  leave  the  main  sequence  first;  thus  it 
will  also  be  the  first  to  reach  its  Roche  limit.  From  the  onset  of 
nuclear  burning  to  the  time  it  fills  its  lobe,  the  system  will  be 
detached.  Both  components  are  in  the  main  sequence,  although  the 


primary  may  be  slighcly  above  for  a  short  while,  from  hydrogen 
depletion  until  it  fills  its  lobe. 

The  primary's  expansion  will  be  checked  at  its  Roche  limit. 

With  further  expansion,  it  will  lose  mass  to  the  secondary.  If 

both  are  still  on  the  main  sequence  at  this  point,  the  transfer 

will  leave  the  primary  the  less  massive  star,  the  new  secondary. 

The  system  is  now  semi-detached  with  a  contact  secondary.  This  is 

called  case  A  mass  transfer,  and  was  depicted  in  Figure  1-8.  If 

the  primary  is  off  the  main  sequence  when  it  fills  its  lobe,  it 

will  lose  its  distended  outer  layers,  which  have  very  low  density. 

Thus  it  loses  much  volume  but  little  mass;  the  mass  ratio  remains 

about  the  same,  but  the  secondary's  Roche  lobe  now  becomes  full. 

After  losing  these  outer  layers,  the  primary  recontracts  to  within 

its  lobe.  This  is  termed  case  B  mass  transfer,  and  again  leads  to 

a  contact  secondary.  Observation  bears  out  this  expectation, 

although  there  is  a  chance  to  observe  a  contact  primary.  This  is 

highly  Improbable,  however,  since  the  time  spent  in  the  contact 

secondary  configuration  is  at  least  two  orders  of  magnitude 

greater  than  its  transient  stay  in  the  contact  primary  con- 
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fiuration.  The  contact  secondary  moves  above  the  main  sequence 
with  the  characteristics  of  a  subgiant. 

The  mass  transfer  leading  to  semi-detached  systems  also  pro¬ 
vides  some  interesting  insights.  For  example,  a  system  with  ini¬ 
tially  highly  unequal  masses  will,  after  mass  transfer,  have  a 
higher,  but  still  low,  mass  ratio.  If  multiple  transfers  occur, 
each  one  will  raise  the  mass  ratio.  This  suggests  that  systems 
with  very  low  mass  ratios  are  young  systems,  although  we  cannot  be 
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certain  that  systems  with  high  (f»i)  mass  ratios  are  old,  since 
they  may  have  formed  in  that  state. 

Theories  concerning  contact  binaries  are  much  more  tenuous. 
Observations  show  that  these  systems  have  minima  of  approximately 
equal  depth,  implying  either  equally  luminous  stars,  a  masking 
layer  of  clrcumstellar  material,  or  both.  Also,  in  all  systems 
studied  for  sufficient  lengths,  the  period  has  decreased  con¬ 
sistently,  which  can  be  explained  by  a  constant  mass  loss  of  about 
10  solar  masses  per  year. 

,  Most  importantly,  contact  binaries  considerably  outnumber  all 
ocher  types.  There  are  two  fundamentally  different  ways  of 
explaining  this  preponderance.  The  first  asserts  that  both  com¬ 
ponents  have  filled  their  Roche  lobes  by  expansion,  have  become 
stable  through  mass  loss  across  an  outer  Lagragian  (points  L2  or 
L3  in  figure  1-7),  and  remain  in  an  evolutionary  trap.  Thus  con¬ 
tact  binaries  would  represent  an  old  configuration.  Observations 
show,  however,  a  marked  absence  of  contact  binaries  in  globular 
clusters  and  other  regions  of  older  star  population.  This  effec¬ 
tively  casts  doubt  on  this  theory  of  contact  formation,  unless 
there  is  some  unspecified  outside  disrupting  force  which  acts  on 
them  in  the  IQ10  year  lifetime  of  old  clusters,  or  unless  the 
regions  of  older  stars  were  uncondusive  to  binary  formation. 

We  can  also  theorize  that  contact  pairs  are  young,  that  they 
form  in  or  nearly  in  their  present  configuration.  This  easily 
explains  their  appearence  in  young  clusters,  and  is  completely 
independent  of  any  future  catastrophes  which  may  cause  their 
absence  from  older  clusters.  However,  their  formation  must  be  an 
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extremely  common  event  to  account  for  their  large  numbers.  If 
true,  this  leads  to  general  predictions  concerning  the  formation 
of  all  stars;  the  contraction  of  extremely  large  pre-stellar 
clouds  about  multiple  mass  centers  should  be  the  rule  rather  than 
the  exception. 

Regardless  of  their  questionable  origin,  they  do  exist  and 
there  are  some  solid  physical  characteristics  of  contact  pairs. 
Since  both  fill  their  Roche  limits,  the  more  massive  must  also  be 
the  largest,  but  the  secondary  can  be  shown  to  be  the  brighter. 

The  existence  of  great  tidal  forces  would  ensure  an  equality  in 
macroscopic  features;  a  system  with  highly  unequal  masses  would 
tear  itself  apart  at  such  close  distance.  Specific  component 
qualities  such  as  spectral  and  luminosity  class  must  await  obser¬ 
vational  placement  of  the  H-R  diagram,  but  we  would  expect  both 
stars  to  be  quite  near  the  main  sequence. 

This,  then,  is  the  theoretical  knowledge  necessary  to 
interpret  the  results  of  this  paper,  the  temperature  and  lumin¬ 
osity  of  individual  components  of  binary  systems,  in  terms  of 
binary  evolution.  The  next  section  will  treat  the  theory  behind 
the  actual  determination  of  these  quantities  from  data  obtained  at 


the  telescope 
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II  -  Experimental  Design 

In  light  of  the  preceding  theoretical  considerations,  one  can 
infer  quite  a  bit  about  stellar  and  binary  evolution  if  one  first 
knows  the  temperatures  and  magnitudes  (i.e.  positions  on  the  H-R 
diagram)  of  the  components  of  binary  systems.  This  therefore  will 
be  the  purpose  of  this  paper:  to  locate  as  many  individual  com¬ 
ponents  of  binary  systems  as  possible  on  the  H-R  diagram  and  to 
investigate  the  consequences  of  this  plot  In  terms  of  binary  evo¬ 
lution.  To  begin,  however,  we  must  see  how  this  information  is 
derived,  first  for  the  case  of  a  single  star,  and  then  for  the 
more  complex  case  of  an  eclipsing  binary. 

As  we  have  already  seen,  the  intensity  of  starlight  at  two 
given  wavelengths  can  give  a  measure  of  stellar  temperature, 
assuming  the  light  received  obeys  the  Planck,  distribution.  If  we 
introduce  a  third  wavelength  (a  second  color  index),  we  will  find 
that  this  reading  does  not  necessarily  fall  on  the  Planck  curve 
found  from  the  first  color  index.  Clearly,  this  means  that  the 
light  received  does  not  conform  neatly  to  a  Planck  distribution; 
the  third  wavelength  tells  how  the  star  differs  from  a  black  body. 

The  second  color  index  is  thus  a  qualitative  description  of  the 
extra-stellar  absorption  and  the  state  of  the  stellar  atmosphere; 
and  extremely  diffuse  atmosphere  (as  in  a  giant)  will  radiate  with 
different  absorption  and  emission  peaks  than  will  a  normal 
atmosphere.  Together,  the  two  color  indices  give  both  the  tem¬ 
perature  and  the  luminosity  class.  This  suffices  to  locate  a  star 
on  the  H-R  diagram,  without  explicitly  determining  its  absolute  magnitude. 
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The  use  of  a  standard,  well  calibrated  set  of  filters  greatly 

facilitates  the  conversion  of  intensity  levels  into  temperature 

25 

and  luminosity  data.  The  Johnson  and  Morgan  UBV  system  was  cho¬ 
sen  for  this  purpose.  This  is  a  set  of  three  filters:  an  ultra¬ 
violet  pass  filter  centered  at  3600A,  a  blue  pass  filter  centered 
at  4250A,  and  a  visible  (yellow-green  pass)  filter  centered  at 
5350A.  There  are  many  other  possible  filter  systems,  but  this  is 
the  most  comprehensively  documented,  and  since  it  is  a  broad  band 
system  (the  three  filters  have  bandwidths  of  3S0A  to  500A),  obser¬ 
vation  times  of  dim  stars  can  be  kept  to  a  minimum,  as  more 
light  passes  through  the  filter.  Tables  for  converting  UBV  data 

into  spectral  and  luminosity  class  (temperature  and  luminosity) 

*  26 
can  be  found  in  Novotny  and  Allen. 

In  order  to  determine  this  information  for  a  binary  system, 

it  must  be  eclipsing,  and  further,  it  must  undergo  totality.  If 

we  observe  the  system  at  totality,  we  will  receive  light  only  from 

the  occulting  star,  and  can  easily  determine  the  color  Indices 

(B-V)  and  (U-B)  for  it.  To  obtain  the  same  results  for  the  other 

star,  which  we  never  see  alone,  we  must  observe  the  system  a 

quarter  phase  out  of  totality,  when  both  stars  radiate  to  us,  and 

then  remove  the  contribution  of  the  known  occulting  star. 

First,  let  us  define  the  magnitude  (at  some  given  wavelength) 

of  both  stars  together,  when  the  system  is  out  of  eclipse: 

mjj  ■  mi&m2  ,  m2  ■  m^j~mi 

where  &  signifies  magnitude  addition  and  ~  signifies  magnitude 
subtraction,  both  of  which  are  dissimilar  to  algebraic  addition 
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and  subtraction,  since  magnitude  follows  a  logarithmic  arithmetic. 

To  find  these  operations,  we  recall  (1-2): 

“^■“1  *  2.51ogU  i/ib)  (1-2) 

this  leads  to: 

*b/*l  -  io^ml“nb^2'5  (2-1) 

A 2  -  *b-‘l  ”  {lO^^b^2-5}^  -  *1  (2-2) 

-|10(m1-mb)/2.5  _  ^  (2-3) 

thus : 

*2/1!  -  { io(“1“nb)/2*5  -  1}  -  10(ml_®2)/2.5  (2-4) 

Taking  logarithms  of  both  sides  and  simplifying  gives : 

m2  -  ml  -  2.51og{l0(“l~1“b)/2*5  -  l}  (2-5) 

which  defines  m^-m!.  Applying  this  equation  to  all  three  wave¬ 
lengths  and  taking  the  proper  differences  can  give  us  the  tem¬ 
perature  and  luminosity  of  both  components  in  any  binary  system 
which  displays  totality. 

We  have  yet,  however,  to  correct  for  the  unique  calibration 
of  our  equipment  in  relation  to  the  Johnson  and  Morgan  standards , 
or  for  the  atmospheric  extinction  which  necessarily  occurs ,  espe¬ 
cially  in  such  a  moist  environment  as  Annapolis.  A  recently 
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published  data  reduction  method  by  Harris  et.  al.  was  selected 
for  this  purpose,  primarily  because  it  solves  for  both  instrumen¬ 
tal  transformation  and  atmospheric  extinction  simultaneously,  and 
in  the  process  takes  advantage  of  the  Naval  Academy's  extensive 
computing  facilities.  The  following  argument  summarizes  this  method. 
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Beginning  by  treating  atmospheric  extinction,  we  define  the 
air  mass,  X,  as  the  distance  through  the  atmosphere  starlight  must 
travel  from  the  star's  position  to  the  observer.  Geometry  shows: 

X*secz( l-.QQ12tan2z)  (2-6) 

secz-(sinLsin6+cosLcos6co8h)“l  ( 2-7 ) 

where  L  is  the  observer’s  latitude,  $  is  the  star's  declination, 
and  h  is  the  star's  hour  angle.  This  latter  quantity  is  a  measure 
of  relative  longitude  of  the  star  as  seen  by  the  observer. 

We  now  assume  that  the  readings  outside  the  atmosphere  are 
related  to  the  measured  readings  by: 

(u-b)m  -  (u-b)0  +  lCubx  ( 2-8) 

The  color  index  (u-b)  is  solely  used  in  this  argument.  Any  other 
would  develop  exactly  the  same.  The  extinction  coefficient 
may  further  be  a  slightly  varying  function  of  color,  which  we  can 
represent  in  a  truncated  power  series: 

kub  “  k"ub  +  ^"ub^u_b)o  (2-9) 

Combining : 

(u-b)m  -  (u-b)0  +  k'ubX  +  k"ub(u-b)0X  (2-10) 

Solving  (2-10)  for  the  reduced  color  index  gives: 

.  -  (u-b)m-k  ub*  (2-11) 

l+k"ubX 

Turning  to  the  Instrumental  transformations,  we  again  write  a 
truncated  power  series,  but  this  time  as  a  coefficient  times  the 


published  color  index: 
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(u-b)0  -  YI  +  Y2<U-B)  +  Y3(^-B)2  (2-12) 

Combining  (2-11)  and  (2-12)  gives: 

(u..b)ln  ZJL 0**  -  yi  +  y2(U'B)  +  Y3(U-B)2  (2-13) 

1  +  Ic  ubX 

After  manipulation  and  a  change  of  constants: 

(u-b)m  -  cl  +  c2X  +  C3(U-B)  +  c4X(U-B)  +  cs(U-B)2  (2-14) 

where:  c1  =  Yi  c2  =  k'ub  +  Yilc"ub 

c3  =  12  =4  =  Y2fc"ub 

c5  =  Y3(1+k">ubx^ 

Our  purpose  now  is  to  find  the  coefficients  c^  so  that  we 
may  transform  the  (u-b)m  into  (U-B).  To  do  this,  we  employ  a  set 
of  standard  photometric  stars  for  which  (U-B)  is  known. 

Next,  we  generalize  this  solution  (2-14)  to  accomodate 
multiple  nights  observing  on  the  same  telescope.  Recalling  the 
definitions  for  the  coefficients  c^,  we  see  that  C3  and  C5  remain 
constant  for  the  same  telescope,  and  c4  can  be  assumed  constant 
over  many  nights,  since  k"ub  is  both  small  and  slow  varying.  We 
rewrite  (2-14): 

Rjk  *  d-b)jit  -  {cut  +  c^Xjfc}  -  {c3  +  c4Xjk}(U-B)j  -  c5(U-B)j2 

where  Rjk  is  the  residuals  tensor,  ideally  zero.  The  J-lndex 
tracks  the  star  out  of  n  total  stars,  and  the  k-lndex  tracks  the 
night  out  of  m  total  nights  observing. 

We  add  a  weighting  factor,  as  a  function  of  the  air  mass 
observed  through,  to  obtain: 
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—  **  2 
s  •  1  1  ^ jltRjk 

k-1  j-1 


(2-15) 


which  we  cry  co  minimize. 

We  now  make  another  change  in  variable  co  write: 


where : 


Rjk  “  (“' 

2m+3 

■b)jk  “  1  8izi jk  (2' 

i-1 

8i  “  cli 

Zijk  ■  ^i^  ( 1< i<m) 

8i  -  cl(i-m) 

zijk  “  ®i-m^^jk  (a+l<i<2m) 

82m+l  ■  c3 

z(2m+l)jk  *  (U-B)j 

82m+2  *  CA 

z(2m+2)jk  -  Xjk(U-B)j 

82m+3  *  c5 

z(2m+3)jk  "  (U-B)j2 

This  allows  us  Co  find  che  minimum  of  S: 


ac  ®  *■ 

H  -0-21  I 


(2-17) 


m  n  2m+3 

°“  I  1  “jk  [C«-b)jk  ~  l  hizijk]zpjk 
k-lj-1  1-1 


(2-18) 


11  a,jk(u_'>) jkzpjk  *  Ill^i® jkzi jkzpjk  (2-19) 

kj  ikj 


For  ease  if  notation,  let  us  here  shift  to  tensor  notation, 
def ining : 


Gu  =  8i 

So  5  *p  "  II“jk(u-b)jkZpjk  (2-20) 

Sio  5  zip  “  jkzijkzpjk 


I 
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Thus  finally. 


S0  -  Z^Gy  (2-21) 


or,  since  Che  coefficients  we  seek  are  concained  in  Gy , 

Gy  -  Zyg^Sg  (  2«  2  2) 

AC  this  point,  since  Che  dimensions  of  these  censors  are  more 
chan  ewice  Che  number  of  nights  observed,  Che  power  of  large  digi¬ 
tal  computers  becomes  attractive.  Appendix  II  shows  the  programs 
used  to  calculate  Gy  from  the  observed  data.  Once  we  have  Gy,  we 
can  re-enter  equation  (2-14)  to  find  (U-B)  for  an  unknown  star 
observed  on  some  night  k. 

The  advantages  of  this  method  primarily  lie  in  the  increased 
ease  and  flexibility  it  gives  the  observations  themselves.  First, 
since  all  standards  observed  affect  all  coefficients,  an  error  in 
any  one  will  be  more  effectively  compensated  by  the  others. 

Further,  since  the  quality  of  the  final  results  depends  on  the 
total  number  of  standards  observed,  the  observer  is  not  forced  to 
study  any  given  standard  at  any  given  time,  allowing  for  capricious 
weather  patterns.  Similarly,  nights  of  poorer  optical  quality  can 
still  be  of  some  use  by  giving  them  a  lower  weight. 

Ue  are  now  armed  with  all  the  theoretical  weaponry  needed  to 
convert  raw  data  into  a  real  measure  of  temperature  and  luminosity 
for  the  individual  components  of  a  binary  system.  Any  modifica¬ 
tion  of  these  techniques  will  arise  through  peculiarities  of  the 
equipment  itself.  Let  us  now  describe  this  equipment  and  its 
effect  on  the  data  taking  process. 
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Clearly,  Che  single  most  important  tool  for  astronomical 
research  is  the  telescope,  used  not  to  magnify,  but  to  gather  as 
much  light  as  possible.  The  Naval  Academy's  telescope  is  a  Ferson 
Optics  1216,  a  sixteen  inch  reflector,  located  on  the  roof  of 
Michelson  Hall.  The  telescope  has  console  and  remote  control  of 
Che  coordinate  slew  motors  (both  slow  and  fast),  the  Cassegrain 
focus,  and  the  observatory  dome.  With  the  photometer  head 
attached  at  the  Cassegrain  focus,  the  telescope  has  a  focal  ratio 
of  14:1. 

The  only  Important  operational  deficiency  of  this  instrument 
was  the  inability  of  the  sidereal  drive  system  (which  should  com¬ 
pensate  for  the  rotation  of  the  Earth)  to  track  a  star  for  periods 
greater  than  a  minute,  especially  when  the  telescope  was  posi¬ 
tioned  far  from  the  vertical  or  when  the  ambient  temperature  was 
either  cold  or  quickly  changing.  Proper  balancing  minimized  this 
problem,  but  it  was  still  everpresent.  Fortunately,  however, 
observational  procedure  could  be  altered  so  that  data  taking 
become  insensitive  to  telescope  drift,  albeit  more  time  consuming. 

The  core  of  the  observational  equipment  was  the  Thom  EMI 
Gerncou  Starlight  1  stellar  photometer.  This  photometer  consists 
of  two  fundamental  parts,  the  photometer  head  and  a  digital 
counter.  The  photometer  head  was  attached  at  the  cassagrain  focus 
of  the  telescope,  and  steadied  against  rotational  and  axial 
displacements  by  a  metal  plate  fixed  fast  to  the  telescope  and 
attached  to  the  head  by  four  long  bolts.  The  head  was  further 
secured  in  a  sleeve  extending  from  the  telescope  by  three  set 


screws 
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The  purpose  of  Che  photometer  head  is  to  detect  light  through 
the  telescope  and  to  convert  it  into  a  useful  electrical  signal. 
Diaphragm  and  filter  turrets,  a  photomultiplier  tube,  and  a  power 
supply  constitute  its  major  components.  The  diaphragm  turret 
allows  one  to  select  one  of  five  aperture  sizes,  ranging  from  .01" 
to  .199".  Each  position  has  two  stops:  one  sets  the  aperture  in 
the  center  of  the  field  of  view,  for  observing  the  star  under 
study,  and  Che  ocher  slightly  offsets  the  aperture  to  give  sky 
background  measurements  while  still  maintaining  the  track  on  the  star. 
This  feature  saved  considerable  amounts  of  time  and  increased  the 
reliability  of  the  results  through  the  added  consistency  it  gave 
the  observations. 

The  filter  turret  operates  in  a  similar  manner  to  allow 
selection  of  a  specific  filter.  The  filters  supplied  in  the  pho¬ 
tometer  head  are  as  follows:  ultra-violet,  1mm.  Shott  UG-1; 
blue,  1mm.  Shott  BG-25;  and  visible,  1mm.  Shott  OG-515.  These  are 
Che  filters  most  compatible  with  the  supplied  photomultiplier 
tube.  The  filter  turret  also  has  two  empty  positions  for  other 
filters  if  desired,  which  were  not  used,  and  an  opaque  filter  to 
protect  the  photomultiplier  from  ambient  light  when  not  in  use. 

The  photomultiplier  Cube  supplied  is  Che  EMI  9924A,  which  is 
ideally  suited  to  UBV  photometry.  Its  spectral  range  is  from 
3100A  to  6700A,  with  peak  response  at  4200A.  This  tube  also  does 
not  need  to  be  cooled  to  give  good  results,  so  chat  the  set  up  of 
the  equipment  was  much  simpler  and  faster.  The  dark  count  of  the 
tube  (the  intrinsic  background  noise)  ranged  from  about  seventy 
counts  per  second  at  0°C  to  ninety  cps  at  20°C. 
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The  counter  receives  Che  output  signal  from  the  photometer 
and  displays  a  number  on  its  L£D  digital  readout.  The  obvious 
advantage  of  a  digital  counter  over  the  usual  strip  chart  recorder 
is  the  higher  precision  achievable  for  low  light  levels.  The 

Q 

maximum  light  count  before  counter  overflow  is  10  ,  which  was 
never  reached  in  any  observation.  The  counter  also  has  several 
sophisticated  capabilities,  such  as  interface  with  other  photo- 
meter  heads  or  with  a  microcomputer,  none  of  which  were  needed  nor 
utilized. 

The  counter  may  be  operated  in  two  modes.  It  starts  on  com¬ 
mand  from  the  observer,  but  may  stop  either  after  a  pre-determined 
time  (.01,  .1,  1,  or  10  seconds),  or  when  the  observer  himself 

a 

stops  the  count  manually.  The  first  mode  was  used  throughout, 
with  a  count  time  of  one  second.  This  time  was  selected  over  a 
longer  one  as  a  compromise  between  poor  telescope  tracking,  as 
previously  discussed,  the  desire  to  keep  the  aperture  size  to  a 
minimum,  usually  .02"  or  .0313",  and  better  statistical  behavior 
of  longer  count  times. 

A  correction  must  also  be  made  for  pulses  too  close  for  the 
electronics  to  resolve  as  distinct.  This  "dead  time”  is  calcu¬ 
lated  by  plotting  log(n)  against  magnitude  for  some  dim  stars  of 
known  magnitude,  where  n  is  the  displayed  count.  This  should 
result  in  a  straight  line.  A  bright  star  is  then  observed,  and 
when  plotted,  will  be  below  this  line.  The  dead  time  is: 

d  -  (N-n)/n2  (2-23) 

where  N  is  the  count  predicted  be  the  plot.  The  corrected  count 
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Figure  2—1*  Schematics  of  the  photometer  head  and  counter. 
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aow  becomes: 


N  -  n(l+dn)  (2-24) 

For  our  Cube,  Che  dead  cime  was  calculated  as  lOOnsec.  We  can  see 
Chat  Chls  Is  a  very  small  correction  for  all  but  the  standards  and 
the  brighter  of  the  eclipsing  binaries. 

Figure  2-1  shows  Che  path  of  a  hypothetical  photon  as  it 
enters  the  photometer  head.  The  photon  first  encounters  a  flip 
prism,  which  is  controlled  by  the  observer  to  either  view  the  star 
or  to  measure  its  intensity.  Notice  that  it  is  before  the  filter 
in  the  optical  path,  so  that  one  can  see  and  track  the  dimmest 
stars  without  light  loss  through  the  filters.  The  prism,  in  the 
view  position,  also  provides  a  safeguard  against  photomultiplier 
damage  due  to  ambient  light. 

Having  this  equipment,  we  now  need  specific  stars  to  observe 

with  it.  For  this  project,  two  sets  of  stars  are  required:  one  set 

of  eclipsing  binary  systems  for  which  temperature-luminosity  data 

are  as  yet  unknown,  and  another  of  standard  stars  with  known  color 

indices,  which  control  the  generation  of  the  reduction  tensors 

discussed  previously.  A  list  of  unknown  eclipsing  systems  was 

2  8 

provided  by  Dr.  Richard  L.  Walker  of  the  U.S.  Naval  Observatory, 
Flagstaff  Station.  His  original  list  consisted  of  some  eighty 
systems  brighter  Chan  twelfth  magnitude,  with  a  good  spread  in 
right  ascension.  The  first  condition  arises  as  a  practical  limit 
of  local  visibility,  a  function  of  both  our  equipment  and  local 
observing  conditions.  The  second  condition  ensures  chat  there 
will  be  something  to  observe  regardless  of  the  time  of  the  year. 
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A  lisc  of  standard  stars  was  compiled  by  the  author  from  che  list 

29 

of  UBV  standards  given  in  Johnson  and  Morgan.  Standards  were 
selected  with  two  criteria:  to  ensure  that  their  co-ordinate 
distribution  followed  that  of  the  program  stars,  minimizing  unde¬ 
tected  positional  variation  in  extinction;  and  to  ensure  a 
spectral  and  luminosity  class  distribution  similar  to  one  expected 
for  components  of  eclipsing  binaries  of  any  type,  minimizing  unde¬ 
tected  variation  in  instrumental  transformation. 

Since  the  beginning  of  the  project,  the  list  of  unknown 
systems  has  undergone  several  revisions  for  a  variety  of  reasons. 
The  most  common  of  these  was  too  long  a  period  for  our  time  frame. 
The  maximum  period  was  arbitrarily  set  at  thirty  days.  To 
illustrate  the  rationale  of  this  choice,  consider  that  not  quite 
one  night  in  seven  was  of  observable  quality.  Unobservability 
can  arise  from  many  causes:  a  cloud  cover,  a  bright  moon,  etc.  It 
is  a  landmark  frustration  of  astronomical  research  that  one  cannot 
control  che  laboratory;  one  can  only  entreat  the  Aesir  and  hope 
for  clear  nights.  If  a  system  eclipses  once  each  month,  then  on 
the  average  we  can  see  only  one  every  seven  months.  Take  into 
account  also  that  in  this  time,  the  Earth  would  precess  through 
more  than  twelve  hours  in  ascension,  and  that  there  is  no  guaran¬ 
tee  that  the  eclipse  will  occur  during  hours  of  local  darkness;  it 
clearly  becomes  more  time  efficient  not  to  bother  with  long  period 
systems,  but  to  concentrate  on  completing  others  with  shorter 
periods. 

To  observe  an  eclipsing  system  in  totality,  other  Information 
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is  also  required.  One  must  know  the  system's  ephemeris,  its  sche¬ 
dule  of  primary  eclipses.  Ephemerides  were  found  in  several  sources, 
primarily  Koch30  and  Rafert.31  If  a  system's  ephemeris  has  not 
been  determined,  however,  we  cannot  use  it  since  the  time  of  tota¬ 
lity  is  unknown.  The  amount  of  time  required  to  calculate  an 
ephemeris  obaervationally ,  especially  with  the  capricious  nature 
of  local  weather  conditions,  would  be  prohibitive  to  completing 
the  initial  purpose  of  the  project. 

One  must  also  be  able  to  locate  the  system  in  order  to 

observe  it.  Co-ordinates  for  all  the  systems  were  found  in 
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Kurkarkin  for  I  January  1900.  These  were  precessed  to  1  January 

1983  for  observational  purposes  and  to  1  January  1950  for 

reference  to  the  Smithsonian  Astrophysical  Observatory  Star 
33 

Atlas  in  order  to  prepare  finder  charts.  The  method  of  pre- 

34 

cession  used  was  the  programmable  calculator  method  of  Shudde. 

Three  systems  were  not  found  on  the  SAO  atlas. 

However,  the  most  limiting  factor  governing  which  systems 
were  observable  was  nothing  so  subtle  as  the  above,  but  rather 
that  the  photometer  system,  which  was  ordered  in  April,  did  not 
arrive  until  mid-January.  Adapting  the  telescope  mountings  to 
compatibility  with  the  photometer  head  and  initial  testing  tied  up 
the  equipment  until  late  January.  Appendix  III  lists  the  observable 
systems  from  Or.  Walker’s  list,  their  co-ordinates,  period,  maximum 
hour  angle  to  allow  an  air  mass  X<1.5  (as  recommended  by 
Harris35),  and  the  number  of  totalities  exhibited  in  local  hours 
of  darkness  within  h-th^j  between  the  date  of  equipment  readiness 
and  the  end  of  April.  One  can  easily  see  that  there  are  several 
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systems  completely  unobservable  and  that  only  about  a  dozen  have  a 
better  than  fifty-fifty  chance  of  showing  just  one  totality,  given 
a  one  in  six  "good  night"  ratio. 

Although  this  severely  limits  the  statistical  nature  of  the 
results,  accuracy  can  still  be  expected  if  observations  are  done 
such  that  each  is  statistically  valid.  To  this  end,  observation 
of  any  event  consisted  of  at  least  ten  one-second  counts  per 
filter,  and  a  standard  star,  preferably  one  near  the  unknown 
system,  was  observed  Immediately  before  and  after.  Also,  somewhat 
paradoxically,  this  forced  the  decision  concerning  the  accep- 
tibility  of  observing  conditions  to  become  more  discriminative; 
more  than  one  observation  has  been  ruined  by  nothing  more  than  a 
small,  stray  cirrus  cloud. 
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III  -  GENERATION  OF  RESULTS 

The  nose  useful  illustration  of  the  method  by  which  values  of 
component  temperature  and  luminosity  are  obtained  would  be  to  take 
a  specific  system  and  trace  the  procedure  step  by  step  to  completion. 
For  this  purpose,  let  us  choose  W  Ursae  Majoris,  the  first  system 
for  which  a  solution  was  obtained.  The  co-ordinates  of  this 
system  are:  9h42m34s  in  right  ascension  and  56°01'49”  in  declina¬ 
tion.  Its  maximum  visual  magnitude  (m*,)  is  46.3,  and  its  magni¬ 
tude  in  totality  (m^)  is  +9.2. 

The  first  step  is  to  calculate  a  schedule  of  total  eclipses 
for  the  system.  Its  ephemerls  is: 

2,434,489.44286  +  .3336384E  (3-1) 

The  first  term  is  a  fixed  time,  an  epoch,  expressed  in  the  astro¬ 
nomical  system  of  Julian  dates,  and  the  second  term  is  the  period, 
in  days,  multiplied  by  an  integral  number  of  eclipses  (E)  to  show 
periodicity. 

The  computer  program  (see  Appendix  II)  written  to  print  out  a 
schedule  of  observable  events  uses  the  first  totality  in  October  1982 
as  a  starting  point.  For  this  ephemerls,  this  is  02:33:50  EST  1 
October.  The  program  then  prints  all  events  for  this  ephemerls 
between  20h00  and  04h00  (local  hours  of  darkness),  assuming  cir¬ 
cular  orbits,  so  that  totality  and  quadrature  are  one-quarter 
period  apart.  The  effect  of  this  assumption  is  quite  valid,  even 
if  the  assumption  itself  is  not,  since  the  stars  are  out  of 
eclipse  for  at  least  a  sixth  of  a  period  before  re-eclipsing. 

Tidal  forces  in  systems  as  close  as  those  to  be  studied  would  easily 
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have  absorbed  any  eccentricities  as  large  as  this. 

There  are  still  two  possible  downfalls  to  the  computer¬ 
generated  event  schedule.  The  first  is  simply  drift  in  the  event 
times  caused  by  round-off  error  in  the  computer.  This  is  easily 
avoided  by  using  double-precision  arithmetic,  accurate  to  twelve 
decimal  places.  The  second  arises  because  the  ephemeris  schedules 
primary  eclipse,  whereas  we  are  interested  only  in  total  eclipses. 
Along  with  his  list,  however.  Or.  Walker  also  indicated  what  kind 
of  eclipse  the  primary  was,  either  annular  or  occulting  (total). 

If  the  primary  eclipse  is  annular,  one  just  adds  half  a  period  to 
the  ephemeris  epoch,  and  uses  this  time  as  a  starting  point  for 
the  program.  If  the  primary  eclipse  is  total,  this  step  is  not 
needed.  The  program  outputs  the  time  (in  EST)  of  an  event  and  its 
corresponding  sidereal  time.  This  is  a  time  reckoned  by  prime 
meridian  passage  of  the  fixed  stars  rather  than  the  sun,  and  is 
entirely  local  (no  system  of  time  zones).  Sidereal  time  allows 
the  calculation  of  hour  angle: 

h  -  ST  -  a  (3-2) 

If  our  sidereal  time  and  a  star's  right  ascension  are  equivalent, 
then  there  is  zero  hour  angle,  and  the  star  is  on  our  meridian. 

We  now  recall  Harris'  previously  discussed  caveat  of  limiting 
the  maximum  air  mass  to  X<1.5.  This  becomes  especially  important 
considering  the  nature  of  this  project;  since  we  have  a  very 
limited  number  of  observations  per  system,  the  accuracy  of  each 
becomes  exceedingly  crucial.  From  equation  (2-6),  the  air  mass  is 
a  function  of  only  the  observer's  latitude,  and  the  star's  declination 
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and  hour  angle.  The  latitude  is  fixed  for  all  observations  at 
38°58',  and  for  a  given  star,  declination  is  also  constant.  Thus 
to  limit  the  air  mass,  we  must  limit  the  hour  angle.  Solving  (2-6) 
first  in  terms  of  cosz  for  h-1.5: 

cos2z( . 9988-1. 5co8z)  *  .0012  (3-3) 


Then  in  terms  of  cosh: 


cosz  -  sinLsinfi 
cos^max  cosLcosfi 


(3-4) 


This  uniquely  determines  a  maximum  hour  angle,  and  events  whose 
hour  angle  is  greater  than  this  are  scratched  from  the  schedule. 

The  next  step  is  to  observe  events  on  our  revised  schedule. 

The  system  is  sighted  by  use  of  the  finder  telescopes  by  comparing 
stellar  field  patterns  to  the  finder  chart.  This  step  was  the 
most  surprisingly  difficult  of  the  entire  observational  procedure, 
since  star  fields  never  quite  seemed  to  agree  exactly  with  those 
published  on  the  SAO  atlas.  Once  the  system  is  located,  it  is 
centered  in  the  smallest  aperture  which  allows  moderately  good 
tracking,  (the  star  remains  in  the  aperture  for  at  least  half  a 
minute).  At  least  ten  readings  are  made  for  each  filter  in  the 
following  sequence:  system  observation  in  V,  background  count  in 
V,  system  and  background  in  B,  then  both  in  U.  Albeit  lengthier 
because  of  increased  setting  changes,  this  cycle  gives  more 
accurate  results  than  other  possibilities,  such  as  taking  all  V 
counts  at  once,  since  time-dependent  errors  (turbulence  or  sparse 
stratus  clouds)  will  be  spread  evenly  over  all  wavelenghts,  mini¬ 
mizing  uncompensated  error  in  any  one.  The  readings  for  the  system 


and  its  background  at  each  wavelength  are  then  averaged,  and  are 
corrected  for  instrumental  dead  time  (2-24).  The  corrected 
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background  count  is  now  subtracted  from  the  corrected  system  count 
to  give  the  contribution  of  just  the  system  itself:  n v,  %,  Oy. 
Color  indices  are  then  computed  via: 

(b-v)  -  -2.51og  Oy/nij  (3-5) 

here  using  (b-v)  as  an  example. 

Before  discussing  the  reduction  of  the  observations  of  the 
program  stars,  let  us  digress  to  the  generation  of  the  reduction 
tensors  used  in  Harris'  method.  Observational  data  of  standard 
stars  is  fed  into  the  computer  (see  Appendix  IV),  which  calculates 
Gy  in  the  sequence  shown  in  Figure  3-1 .  It  is  Important  to  note 
here  the  precision  of  the  inversion  of  Zy0;  the  accuracy  of  Gy , 
and  thus  the  validity  of  the  entire  set  of  results  is  at  the  mercy 
of  this  single  operation. 

With  Gy ,  we  are  ready  to  find  the  actual  color  indices  of  the 
system.  Equation  (3-6)  is  the  governing  relation  for  this 
process : 

0  -  (b-v)  -  glk  -  g2kXk  -  (g3  +  g4Xk)(B-V)  -  g5(B-V)2  (3-6) 

This  is  exactly  equation  (2-16),  but  with  the  residuals  tensor, 

Rjk,  set  equal  to  zero.  Taking  the  proper  elements  of  Gy ,  it  is  a 
simple  matter  to  calculate  the  actual  color  index,  (B-V). 

To  find  the  proper  elements  of  Gy ,  we  recall  its  theoretical 
derivation  (2-20).  Gy  has  2sH-3  elements,  where  m  is  the  total 
number  of  nights  observed.  The  night  of  the  specific  observation 
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♦  OLD  YLEYN 

«aaa«««fl  <~  yleyn  ; 

*ituN  6A5IC7 


I N  P  U  T  #  Or  NIGHTS.  # 
INPUT  1  r  UK  <S-V)f  2 


pas  swore 

OF  STARS'  1 0 .  i Z 
Fuit  iu'B;  >  3  Fuft 


w — v  ?  Z 


*OLD  TEN50P 

aaaaaaaa  < —  “tensop"  password 

*RUN 

INPUT  #  uF  NIGHTS?  10 

DGEST  THOU  NISH  TO  SEE  Z  AND  5?  NAY 


DQE57  THOU  NISH  TO  SEE  Zi-li?  NAY 


DuEST  THOU  NISH  THE  ACCURACY  OF  Zi-i)  CONFIRMED”  vERI 
Z .  77703E-Z  IS  THE  SUM  OF  TuTAl.  ABSOLUTE  DEvIATIui\  OF 
FROM  THE  IDENTITY  TENSOR. 


D0E5T  THOU  NISH  TO  SEE  Z*Z(-i>?  NAY 
*•**  FIAT  0  •*** 


0 . 553246 
0. i 56366 
0. 753635 
-3. 17451 
-1 . 0073i 
1.07757 
3.7S274 
1 .45077 
3 . 30S35 
1.2207S 
0.651254 
1 . 46005 
0 . 533666 
4 . 71305 
2 . 5665 
0.551152 
-1.53425 
0.167145 
-2.00736 
0.312505 
1 . 5225 
-0 . 520i44 


-7. 612Z5E-2 


i_  < 

Z-*Z «.  -  i  i 


Figure  3-1.  Generation  of  the  G  tensor 


£  ^ 

secs  k;  thus  gilt  is  Che  k  element,  g^  is  the  (m+k)  element, 
and  g3,  g4,  and  g5  are  the  last  three  elements.  As  an  aside,  we 
can  see  that  truncating  our  power  series  (2-9), (2-12)  was  indeed 
justified  by  the  small  co-efficient  g5« 

Inserting  the  data  for  U  Ursae  Majoris  into  (3-6)  gives  the 
following  results : 


(B-V)b  -  .590  (B-V)i  -  .584 

(U-B)b  -  .108  (U-B)i  -  .175 


(3-7) 


The  tables  in  Novotny37  use  (U-V)  as  an  entering  index,  which  is 
easily  calculated: 


(U-V)  -  (U-B)  +  (B-V)  (3-8) 

We  can  immediately  find  the  spectral  and  luminosity  class  of 
the  occulting  star,  which  will  locate  it  on  the  ft-R  diagram,  by 
using  the  forementioned  tables. 

To  determine  this  information  for  the  second  star  requires 
more  work.  We  must  subtract  the  contribution  of  the  occulting 
star  from  the  measurements  made  outside  eclipse,  but  we  must  do 
this  via  (2-5): 

m2  -  mi  -  2.51og{  io(mi_1"b)/2*5  -  l}  (2-5) 

Since  a  color  index  is  merely  the  difference  in  magnitude  at  two 
different  wavelengths,  we  adapt  this  to: 

(b-v)2  -  bx  -  2.51og{l0(bl-bb)/2,5  -  1} 

-  vL  +  2.51og{l0(vi"Vb)/,2‘5  -  1}  (3-9) 


(b-v>2  *  (b-v)i  +  2.51og 


(3-10) 


-1Q<vrvb>/2-5  -  n 

^Q(bj^-bb)/2.5  _  ^  | 

Equation  (3-10)  is  in  the  proper  form,  but  to  use  it,  we  need  a 
system  of  maginitudes.  We  could  use: 

v  -  -2.51ognv  (3-11) 

where  Oy  is  the  count  in  the  visible  band.  This,  however,  in  no 
way  accounts  for  variations  in  observing  conditions  among  nights, 
or  any  other  of  the  possible  effects  on  any  single  observation. 

The  form  of  (3-11)  is  still  valid;  but  we  oust  first  normalize  our 
counts . 

First,  we  taite  Qyb,  the  visible  count  of  the  system  out  of 
eclipse.  We  then  assign  n^  and  n^  by  the  relation: 

Uyb  -  10(b-V)b/2*5nbb  (3-12) 

where  (b-v)b  is  the  result  of  applying  (3-6),  given  the  indices  (3-7). 

To  assign  counts  for  the  occulting  star  alone,  we  first  find 
the  brightness  of  it  in  terms  of  the  brightness  of  both: 

nVl/nVb  -  iVl/£Vb  -  IO(n,l_Bb)/2*5  (3-13) 

We  would  expect,  that  if  the  occulting  star  is  one-half  as  bright  as 
both  stars  together,  it  would  have  half  as  many  counts.  Visual 
magnitudes  for  the  systems  were  found  primarily  in  Koch.  Once  we 
define  Oy^,  the  others  are  found  by  again  using  (3-12).  For  W 
Ursae  Majorls,  these  counts  are  listed  below: 
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nvb  -  «49 

nVl  -  2282 

nbb  *  2526 

nbl  -  1333 

(3-14) 

nUfe  -  2287 

nUl  -  1134 

To  obtain  a  consistent  magnitude  system,  we  take  the  logarithms  of 
these  normalized  counts,  but  we  can  dispense  with  the  factor  of 
2.5,  since  it  also  appears  in  the  denominator  of  the  exponent  of 
(3-10).  For  W  Ursae  Majoris,  (3-10)  gives: 

(B-V)2  -  .597  (U-V)2  -  .657  (3-15) 

Before  entering  the  tables  with  the  values  in  (3-7)  and 
(3-15),  let  us  comment  on  the  statistical  inaccuracy  of  Harris' 
method,  introduced  by  setting  Rj^  to  zero.  This  is  the  case  for 
the  aggregate  of  standard  star  data,  but  is  not  necessarily  true 
for  an  individual  star.  Using  standard  stars  as  though  they  were 
unknowns,  (3-6)  gives  results  that  are  slightly  off  the  published 
color  Indices.  However,  we  can  use  (B-V)  and  (U-V)  independently 
to  obtain  temperature  and  luminosity,  and  then  assign  these  two 
pairs  of  values  as  the  endpoints  of  an  uncertainty  bar.  In  this 
way,  we  can  have  a  graphical  measure  of  the  precision  of  the 
results  without  having  to  subjectively  assign  error  bars  based  on 
the  behavior  of  the  standard  stars. 

For  U  Ursae  Majoris,  the  tables  tell  us  that  both  stars  are 
F9V  stars.  The  secondary's  locus  on  the  tt-R  diagram  runs  from 
(-3.776,0.3)  to  (-3.775,0.22)  in  co-ordinates  of  (-log  Teff,log  L/Lp). 
The  locus  of  the  primary  is  (-3.778,0.065)  to  (-3.789,0.10).  We 
can  see  that  these  results  are  at  least  reasonable;  the  system's 
period  suggests  that  it  is  a  contact  binary,  which  is  supported  by 
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Che  poscioas  of  ics  components  on  Che  tt-R  diagram,  recalling  che 
expectations  for  a  contact  pair  from  Che  first  section. 

The  alert  reader  will  have  noticed  that  we  have  not  mentioned 
the  light  varying  effects  also  discussed  in  the  first  section.  A 
quick,  examination  of  these  will  show  why  they  can  be  ignored. 

Host  importantly,  since  we  observe  the  system  only  at  two 
instants,  totality  and  quadrature,  their  time  dependence  becomes 
irrelevant.  Distortion  now  only  makes  the  occulting  star  radiate 
to  us  through  a  smaller  cross-sectional  area,  and  makes  the  stars 
at  quadrature  radiate  through  a  greater  area.  However,  the 
resulting  differences  in  instantaneous  brightness  is  already  taken 
into  account  in  the  measured  visual  magnitude  of  the  system  in  and 
out  of  eclipse,  as  found  in  Koch. 

We  can  dismiss  much  of  the  complications  of  darkening  mech¬ 
anisms  by  a  similar  argument.  Since  we  observe  only  at  totality 
and  quadrature,  we  never  see  one  star  crossing  the  limb  of  the 
other,  where  the  darkening  mechanisms  have  their  greatest  effect 
on  the  light  curve.  Specifically,  in  the  quasi-static  limit,  limb 
darkening  would  make  both  scars  appear  slightly  cooler  than  they 
really  are,  both  in  and  out  of  eclipse.  Gravity  darkening  will 
make  the  occulting  star  cooler  and  dimmer  yet  in  totality,  but 
will  make  both  stars  appear  relatively  brighter  out  of  eclipse. 
Thus  the  overall  effect  would  be  a  solution  for  the  occulting  scar 
slightly  too  cool  and  dim,  and  since  limb  and  gravity  darkening 
offset  each  other  at  quadrature,  the  solution  for  the  other  comp¬ 
onent  would  be  slightly  too  warm.  These,  however,  are  insignifl- 
gant  compared  to  the  uncertainties  discussed  above. 
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The  irradiation  effect  is  completely  irrelevant  at  totality, 

since  all  reflected  or  re-emitted  light  travels  away  from  us.  Out 

of  eclipse,  there  will  be  an  irradiation  eifect,  but  we  will 

asstn&e  that  it  is  also  a  negligible  factor,  as  Dr.  Walker  recom- 
3  8 

mends,  in  order  to  process  as  much  data  as  possible.  The  slight 
increase  in  accuracy  would  not  merit  the  greatly  Increased  time 
required  to  reduce  this  effect,  which,  like  the  darkening  mech¬ 
anisms,  would  lie  well  within  the  uncertainty  bars. 

We  have  now  gained  all  the  information  which  we  set  out  to 
determine  for  this  system;  we  know  the  spectral  class  of  each  com¬ 
ponent,  which  gives  a  measure  of  temperature.  More  importantly, 
the  position  of  each  component  on  the  tt-R  diagram,  which  gives 
clues  to  the  evolutionary  process  in  binaries.  The  procedure 
outlined • above  can  be  used  to  obtain  the  same  information  for  any 
eclipsing  binary  system  exhibiting  totality.  The  next  section 
will  highlight  the  important  points  of  this  procedure  for  the 
systems  observed,  display  the  results  graphically  on  H-R 
diagrams,  and  discuss  the  slgnifigance  of  these  plots. 
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IV  -  DATA  AND  INTERPRETATION 

The  data  taken  for  any  event  consists  of  about  sixty  inten¬ 
sity  measurements  (ten  readings  over  three  filters  at  two  settings 
each)  and  an  average  time  of  observation.  The  value  of  presenting 
this  information  in  this  format  is  questionable  at  best,  as  it 
would  serve  only  to  engulf  the  reader  in  unnecessary  detail.  One 
would  become  detached  from  the  more  important  process  of 
generating  the  final  results. 

To  keep  this  sense  of  flow  in  the  data  analysis,  values  used 
in  key  stages  of  the  result  generation  algorithm  discussed  in  the 
preceding  section  will  be  given  for  all  systems  observed.  The 
data  page  for  a  system  will  therefore  begin  with  basic  information 
about  the  system,  coordinates  and  the  ephemerls  used. 

Observational  data,  after  being  corrected  for  instrumental  dead 
time,  will  then  be  given,  followed  immediately  by  the  resultant 
reduced  color  Indices.  Aa  an  Intermediate  step  in  the  calculation 
of  the  indices  of  the  occulted  star,  the  normalized  counts  will 
also  be  given.  Finally,  once  we  know  the  indices  for  both  com¬ 
ponents,  their  spectral  and  luminosity  classes  will  be  given. 

An  H-R  diagram  showing  the  positions  of  the  components  will  imme¬ 
diately  follow  each  system's  data  page. 

The  main  sequence  on  each  of  these  diagrams  is  the  locus  of 
points  for  the  luminosity  class  V  stars  from  the  tables  in 
Novotny.  The  three  diagrams  immediately  following  illustrate 
various  features  of  this  main  sequence.  The  first  shows  the  spe¬ 
cific  points  from  the  table  which  define  the  main  sequence.  Because 
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we  are  dealing  with  discrete  points,  the  band  is  not  smooth.  This 
will  become  more  evident  later.  The  second  plots  the  sun,  pri¬ 


marily  to  be  used  as  a  reference  for  future  diagrams.  Already 
familiar  with  the  characteristics  of  the  sun,  we  can  gain  a  idea 
of  the  physical  states  of  the  components  in  an  unkown  system  by 
comparing  their  locations  to  the  sun's.  The  third  H-R  diagram  is 
an  expansion  of  the  second  about  the  area  near  the  sun.  This 
clearly  shows  the  jaggedness  of  our  main  sequence.  This  plot 
points  out  another  of  its  shortcomings.  The  sun  we  know  to  be  on 
the  main  sequence,  but  on  this  plot,  it  is  not.  This  is  because 
we  have  represented  a  band  by  a  line.  The  reader  whose  aesthetic 
balance  is  upset  by  these  inaccuracies  may  mentally  broaden  and 
smooth  the  line  representing  the  main  sequence  on  the  diagrams. 
However,  since  the  same  line  is  used  in  all  plots,  this  is  not  a 
critical  point. 
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TU  CAMELOPARD I 

a  *  5h53m26s  <5  -  59O53'09" 

Ephem:  2432633.6555  +  2.933229E 
Primary  Eclipse  -  Annular 


Totality:  4  Apr 
Out:  2  Mar 


X  V  B 

1.211  37742  34285 

1.077  48822  42313 


(b-v)  (u-b) 

6713  .10  1.77 

8777  .16  1.71 


Reduced  Color  Indices: 

(B-V)h  *  -073  (B-V)!  -  -.004 

(U-V)b  -  .059  ( U-V ) !  -  .022 


mb  -  5.5 


mtot  “  5.6 


Normalized  Counts: 


New  Color  Indices: 


Vb  -  48822 
8b  -  45647 
Ub  ■  46240 


(B-V) 2  -  1.544 


Vi  *  44526 
w‘  44690 
UL  -  43633 


( U-V) 2  *  2.14 


Spectral  Classes: 


Primary  -  A1V 
Secondary  -  K4IV 


CAME LOP* 
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TW  DRACONIS 

a  -  15h33m36s  5  -  63°57'49" 

Ephem:  33888.452  +  2.S0687E 
Primary  Eclipse  -  Occulting 

X  V 

Totality:  21  Apr  1.027  5201 
Out:  16  Apr  1.214  5216 

Reduced  Color  Indices: 

( B-V ) b  »  .193  ( B-V) [  -  1.164 

(U-v)b  ■  .267  (U-V)i  *  2.332 


rnj,  “  7.5  mtot  ”  9. 8 


Normalized  Counts: 

Vb  -  5216  Vy  -  522 
Bb  ■  4367  »  178 
Ub  -  4079  '  Ui  -  61 


New  Color  Indices:  (B-V) 2  ■  *122 


Spectral  Classes:  Primary  -  A3V 
Secondary  -  K2IV 


B  £  (b-v)  (u-b) 

2048  143  1.012  2.89 

4369  845  .192  1.784 


(U-V) 2  -  *166 


IAGRAMt  TU  DRACONIS 
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log  Toff 


DIAGRAM*  COMPOSITE  SEMI-DETACHED 
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W  URSAE  MAJORIS 

a  -  09h42m343  6  -  56°01'49" 

Ephem:  2434489.44286  +  .3336384E 
Primary  Eclipse  »  Occulting 


X 

V 

B 

U 

(b-v) 

(u-b) 

Totality: 

2  Mar 

1.056 

4070 

2450 

428 

.551 

1.894 

15  Mar 

1.111 

3458 

2337 

386 

.425 

1.955 

Out: 

2  Mar 

1.070 

4244 

2516 

484 

.568 

1.790 

2  Mar 

1.181 

4459 

2918 

514 

.575 

1.885 

Reduced  Color  Indices : 

(B-V)b  -  .590  (B-V)l  -  .584 

(U-V)b  -  .698  (U-V)L  -  .759 


mb  -  8.5  ®tot  ”9.2 


Normalized  Counts: 

Vb  -  4349  Vi  -  2282 
Bb  -  2526  Bi  -  1333 
Ub  -  2287  Ui  -  1134 


New  Color  Indices:  (B-V>2  ■  .597  (U-V)2  ■  .657 


Specral  Classes:  Primary  -  GOV 

Secondary  -  F9V 
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V566  OPHIUCHI 


a  -  I7h56m02s 

5  -  4°  59 '  16" 

Ephem :  2435245.544 

+  .4096410  IE 

primary  Eclipse  ■  Annular 

Tocalicy:  16  APr 

Out:  20  Apr 

21  Apr 

X  v 

1.263  5526 

1.250  4442 

1.219  4751 

B 

3921 

3066 

2031 

U 

816 

666 

595 

Reduced  Color  Indices: 

(B-V)b  »  .461 
(U-V)b  -  .488 

(E-V)i  -  .449 
(U-V)i  -  -523 

mb  ■  7.5 

mtot 

-  7.9 

Normalized  Counts : 

Vb  -  4597 
Bb  -  3007 
Ub  -  2933 

VL  - 
Bi  - 
Ul  > 

3180 

2103 

1964 

New  Color  indices: 

(B-V) 2  * 

.490 

(U-V) 2  - 

Spectral  Classes: 

Primary  - 
Secondary  - 

F6V 

F5.5V 

(b-v) 

.373 

.403 

.468 


.413 


l 


I 

I 


(u-b) 

1.784 

1.658 

1.768 


AGRAM 


DIAGRAM I 
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AK  HERCULIS 


a  -  17hl3ml2s  <5  -  16°22' 15" 
Ephem:  2436757.6601  +  .42152502E 
Primary  Eclipse  »  Annular 


X  V  B  V  (b-v)  (u-b) 

Totality:  20  Apr  1.149  2007  1283  256  .486  1.75 

Out:  16  Apr  1.417  2022  1458  279  .413  1.795 


Reduced  Color  Indices : 

( B-V ) b  -  .459  ( B-V) j  -  .536 

(U-V)b  -  .544  (.U-V)L  -  .607 


IAGRAMl  A K  HERCULIS 


log  Toff 


H-R  DIAGRAM*  AK  HERCULIS 


log  T«ff 


H-R  DIAGRAM*  COMPOSITE  CONTACT 


log  Toff 
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AG  VIRGINIS 

a  -  12h00mlls  5  -  13°06'10" 

Ephem:  2434086.4176  +  .64264907E 
Primary  Eclipse  -  Annular 


X 

V 

B 

U 

(b-v) 

(u-b) 

Totality: 

16 

Apr 

1.094 

1502 

1320 

261 

.14 

1.76 

Out: 

15 

Mar 

1.124 

2227 

1871 

329 

.189 

1.89 

23 

Mar 

1.224 

2626 

2423 

525 

.087 

1.66 

Reduced  Color  Indices: 

(B-V)b  -  .132  (B-V)i  -  .068 

(U-V)b  -  .193  (U-V)i  -  .127 


mb  -  9.0  mt0C  -  9.4 


Normalized  Counts: 


Vb  -  2427 

Vl 

-  1680 

Bb  -  2149 

Bl 

»  1578 

Ub  -  2032 

«1 

-  1495 

New  Color  Indices:  (B-V) 2  *  .295 

(U-V)2 

Spectral  Classes:  Primary  -  A2V 
Secondary  -  A9V 


IAGRAM I  AG  UIRGINIS 


log  Tmft 
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a  CORONA  BOREALIS 
a  -  15h33m58s  6  -  26O46'20" 

Ephem:  2423163.770  +  I7.359907E 
Primary  Eclipse  ■  Annular 


X  V  B  U  (b-v)  (u-b) 

Totality:  20  Apr  1.076  759610  738445  165824  .031  1.622 


Reduced  Color  Indices:  (B-V)^  -  -.108 

(U-V)L  -  -.169 


Spectral  Class:  Primary 


B9V 
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AM  LEONIS 

a  -  UhOlmlSs  6  -  9°59'21" 

Ephem:  2435593.7168  +  .36579912E 
Primary  Eclipse  -  Occulting 

X  V  1  £  (b-v) 

Totality:  12  Apr  1.201  864  553  81  -484 

Reduced  Color  Indices:  (B-V)i  -  .485 

(U-V)i  -  .633 

Spectral  Class:  Secondary  -  F7V 


(u-b) 

1.983 


H-R  DIAGRAM l  a  CrB  AND  AMLEO 
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TZ  BOOTES 

a  -  15h07m32s  6  -  40°02'01" 

Ephem:  2432701.428  +  .297 16E 


Out: 


X  V 

23  Mar  1.003  2846 


JJ 

564 


U  (b-v) 

96“  1.76 


Reduced  Color  Indices:  (B-V)b  *  2.067 

(U-V)b  -  2.325 


RS  CANES  VENATICI 

a  -  13a09m50s  6  -  36°01'29" 

Ephem:  2433016.819  +  4.79781E 


V  _B  _U  (b-v) 

35?9  2071  387  .588 

(B-V)b  -  . 

(u-V)b  ■  • 


X 

Out:  12  Apr  1.012 


(u-b) 

1.92 


(u-b) 

1.821 


Reduced  Color  Indices: 


632 

680 
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The  Hr-R  diagrams  of  TU  Camelopard!  and  TW  Draconls  strongly 
suggest  that  they  are  semi-detached  systems.  Both  primaries  are 
relatively  large  main  sequence  scars  (spectral  classes  A1V  and 
A3V),  and  both  secondaries  are  above  the  main  sequence,  as  we 
would  expect  for  a  contact  secondary  with  a  diffuse  atmosphere. 

Both  periods  are  about  2.3  days,  also  what  we  would  expect  for  a 
semi-detached  system. 

Similarly,  the  diagrams  of  W  Ursae  Majorls,  V566  Ophluchi, 
and  AK  Herculls  suggest  that  these  are  all  contact  systems.  The 
large  scale  diagram  shows  Chat  both  components  In  any  system  are 
grouped  around  the  main  sequence.  The  small  scale  diagram  shows 
that  the  secondary  Is  Che  brighter  in  all  three  cases,  although 
both  components  have  about  the  same  temperature.  all  of  the  com¬ 
ponents,  whether  primary  or  secondary,  fall  within  a  very  limited 
spectral  class  range,  from  FIV  to  GOV,  with  most  being  lace  (F6+F9). 
This  uniformity  among  the  components  is  in  agreement  with  theore¬ 
tical  predictions  stemming  from  tidal  forces  as  discussed  in  the 
first  section.  The  periods  of  these  systems,  0.33  to  0.46  days, 
also  suggest  contact  binary  status. 

The  system  AG  Virginia  presents  a  unique  problem.  Its 
period,  .62  days,  is  within  the  expected  range  of  periods  for  con¬ 
tact  systems.  However,  the  spectral  classes  of  its  components, 

A2V  and  A9V,  are  not  what  we  would  expect  in  a  contact  system,  and 
its  primary  is  brighter,  also  in  contradiction  of  our  expec¬ 
tations.  The  components  are  also  of  the  wrong  spectral  classes  to 
be  a  semi-detached  system,  as  their  positions  on  the  H-R  diagram, 
at  a  quick  glance,  may  suggest. 
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Thus  we  are  lefc  with  Che  possibility  of  a  detached  system. 
There  is  nothing  in  the  photometric  elements  of  the  results  which 
would  disallow  this;  we  expect  the  primary  to  oe  hotter  and  more 
luminous  in  a  detached  pair,  and  such  a  system  has  no  real 
restrictions  on  component  spectral  class.  Unfortunately,  the 
period  is  very  short  for  a  detached  pair.  One  method  for 
decreasing  the  period  while  maintaining  separation  is  to  Increase 
orbital  speed.  If,  however,  the  stars  did  revolve  about  each 
other  at  high  speeds,  their  atmospheres  would  become  distended,  as 
though  the  stars  were  leaving  wakes.  This  would  cause  both  com¬ 
ponents  to  move  somewhat  above  the  main  sequence,  which  does  not 
occur  in  the  primary. 

A  more  intriguing  possibility,  although  one  requiring  further 
study  far  beyond  the  time  frame  of  this  project,  is  that  this 
system  is  in  a  state  of  change.  Suppose  the  system  has  recently 
undergone  a  case  8  mass  transfer.  The  primary  has  just  settled 
back  to  the  main  sequence  (a  rapid  process),  while  the  secondary 
is  starting  to  move  off,  due  to  its  newly  acquired  atmosphere. 
Recalling  Figure  1-5,  the  secondary  will  continue  to  move  off  into 
the  giant  or  sub-giant  region,  and  the  system  would  become  semi¬ 
detached.  Further  study  may  bear  out  the  validity  of  this 
hypothesis,  but  even  though  greatly  accelerated  due  to  the  artifi¬ 
cial  prodding  of  mass  transfer,  the  transit  to  the  giant  region 
will  be  measured  in  millennia. 

One  may  similarly  hypothesize  that  the  secondary  is  in  transit, 
not  away  from,  but  to  the  main  sequence.  This,  however,  can 
explain  fewer  observed  features  of  the  system.  Most  importantly, 
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a  pre-main  sequence  scar  at  che  secondary's  poslcion  would  noc 
have  an  overly  diffuse  atmosphere,  so  chac  che  problem  of  che 
shore  period  scill  remains.  Also,  che  sear  could  cake  up  to  10T 
years  Co  seccle  on  che  main  sequence,  so  chac  we  would  noc  have 
che  opporcunicy  .u  verify  if  Chis  is  indeed  occurring  for  quice  a 
while.  The  hypochesis  of  acceleraCed  secondary  eranslc  away  from 
che  main  sequence  following  case  B  mass  cransfer,  leading  co  a 
semi-detached  configuration,  while  noc  without  fault,  can  best 
explain  boch  che  physical  and  photometric  features  of  che  system, 
and  ic  is  also  che  easiest  Co  observatlonally  investigate. 

The  H-R  diagram  for  ex  Corona  Borealis  and  AH  Leonls  shows 
just  one  scar  from  each  system  since  these  were  observed  only  at 
totality.  We  can  predict,  however,  based  on  chelr  periods  and  che 
location  of  one  of  their  components,  to  what  kind  of  system  each 
belongs.  Wich  a  B9V  primary  and  a  long  period,  a  Corona  Borealis 
is  either  a  detiched  or  semi-detached  system.  AM  Leonls,  on  che 
other  hand,  with  an  F7V  secondary  and  a  short  period,  is  in  all 
probability  a  contact  system.  No  information  useful  on  an  H-R 
diagram  could  be  determined  for  TZ  Bootes  and  RS  Canes  Venaticl, 
since  these  systems  were  observed  only  out  of  eclipse. 
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V  -  CONCLUSIONS 


The  individual  temperatures  and  luminosities  of  both  com¬ 
ponents  were  determined  for  six  binary  systems,  as  summarized  below: 


SYSTEM 

PRIMARY 

SECONDARY 

SYSTEM  TYPE 

TU  Cam 

A1V 

K4IV 

Semi-Detached 

TW  Dra 

A3V 

K2IV 

Semi-Detached 

AK  Her 

P8V 

F1V 

Contact 

V566  Oph 

P6V 

F5.5V 

Contact 

W  UMa 

GOV 

F9V 

Contact 

AG  Vir 

A2V 

A9V 

- 

AG  Virginia  does  not  fit  neatly  into  a  particular  system 
type.  Its  photometric  elements  suggest  a  detached  pair,  but  its 
period  seems  too  short  for  such  a  configuration.  It  is  hypothe¬ 
sized  that  the .secondary  is  in  transit  from  the  main  sequence 
following  case  B  mass  transfer,  and  thus  that  the  system  is  in  the 
process  of  becoming  semi-detached.  Additional  study  is  needed  on 
a  long  term  time  scale  in  order  to  substantiate  this  hypothesis. 

Solutions  for  the  occulting  star  were  obtained  for  two  addi¬ 
tional  systems: 

SYSTEM  PRIMARY  SECONDARY 

a  CrB  B9V  - 

AM  Leo  -  F7V 

The  systems  TZ  Bootes  and  RS  Canes  Venatici  were  observed  only  out 
of  eclipse,  so  that  no  definitive  solution  for  either  of  their 
components  could  be  reached.  Complete  solutions  for  these  last 
four  systems  may  be  obtained,  however,  with  an  observation  of  the 
proper  event. 

The  reduction  method  of  Harris,  as  programmed  onto  the  Naval 
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Academy's  computing  facilities,  proved  itself  to  be  quite  effec¬ 
tive  in  transforming  raw  observational  data.  The  shortness  of  the 
uncertainty  bars  on  the  Hertzsprung-Sussell  diagrams  attests  to 
the  precision  of  the  method,  and  the  ease  with  which  solutions 
were  reached,  after  the  computer  programs  had  been  written,  cannot 
be  overstated.  But  more  importantly,  besides  proving  to  be  pre¬ 
cise  and  fast,  the  reduction  seems  accurate,  in  that  it  leads  to 
final  results  consistent  with  theoretical  expectations.  These 
factors  make  this  method  highly  recommendable  to  any  photometric 
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APPENDIX  I 


NUCLEAR  FUEL  CYCLES 

The  three  most  important  fusion  cycles  to  stellar  energy 
generation  are  the  proton-proton  chain,  the  carbon-nitrogen-oxygen 
cycle,  and  the  triple-alpha  chain.  Each  will  be  analyzed,  con¬ 
centrating  on  energy  generated,  fractional  mass  loss  per  cycle, 
and  energy  loss  through  neutrinos. 


I  Proton-Proton  Chain  (T<15X106K) 
a)  2X:  H1  +  H1  ♦  D2  +  e+  +  v 
0)  2X:  E?  +  H1  >  He3  +  y 

y)  He3  +  He3  He4  +  H1  +  H1 


Ee  +  Ey  -  1.4421  MeV 
Ey  -  5.4935  MeV 
E  -  12.8599  MeV 


NET:  4H1  ♦  He4  +  E 


E  -  26.7312  MeV 


This  net  energy  release  is  only  0.7119Z  of  the  original  mass 

of  the  four  hydrogens.  From  this,  we  can  assume  that  the  mass 

during  the  main  sequence  is  constant.  We  can  get  a  rough  maximum 

for  the  neutrino  loss  by  allowing  the  positron  to  rapidly  annihilate 

with  an  electron,  releasing  two  .511MeV  y ' s «  leaving  .4201MeV  to 

the  neutrino.  The  average  neutrino  loss,  however,  is  only  about 
39 

. 26MeV,  or  about  1.945Z  of  the  energy  produced. 


i 


(15*106<T<1G8) 


II  Carbon-Nitrogen-Oxygen  Cycle 

а)  C12  +  H1  ♦  N13  +  y  Ey  -  1.9436MeV 

8)  N13  ♦  C13  +  e+  +  v  Ee  +  Ey  -  2.2205MeV 

Y)  C13  +  H1  ♦  N14  +  y  Ey  -  7.5499MeV 

б)  N14  +  H1  ♦  015  +  y  Ey  -  7,2922MeV 

e)  O15  ♦  N15  +  e+  +  v  Ee  +  Ey  -  2.7606MeV 

Z)  N15  +  H1  ♦  C12  +  He4  E  -  4.9646MeV 

NET:  4Hl  -*>  He4  +  E  Et  -  26.7312MeV 

12 

The  C  acts  as  a  catalyst  allowing  reaction  (a)  to  go. 

Again  the  energy  released  Is  only  .7119Z  of  the  Initial  mass, 

giving  constant  mass  during  this  fusion  process.  Rough  maxima  for 

the  neutrino  loss  are  1.1985MeV  for  reaction  (B)  and  1.7386MeV  for 

reaction  (e).  More  sophisticated  treatment  gives  average  neutrino 

40 

loss  as  .7105MeV  and  1.0006MeV  respectively.  The  total  energy 
loss  to  neutrinos  is  therefore  l.7lllMeV,  or  6.4011Z  of  the  total 
energy.  This  loss  is  clearly  no  longer  negligible,  and  will 
become  an  Important  factor  in  later  evolution,  as  discussed  in 
section  I. 

III  Triple  Alpha  Chain  (T>108K) 

a)  He4  +  He4  ♦  Be8  +  y  Ey  -  -.094641MeV 

B)  He4  +  Be8  ♦  C12  +  y  Ey  -  7.3683MeV 

NET:  3He4  ♦  C12  E  Et  -  7.2736MeV 

8  -l  $ 

Since  de  has  an  extremely  short  half-life,  2X10  s,  there 
must  be  a  very  high  density  of  He  for  reaction  (B)  to  go.  A 
helium  nucleus  whose  kinetic  energy  is  solely  derived  from  the  temperature 
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3 

can  travel  only  about  1.289A  in  one  Be  half-life.  The  energy 
released  amounts  to  only  .065032  of  the  mass  of  the  three  helium. 
Additional  reactions  further  convert  some  of  the  carbon  into  oxy¬ 
gen  to  form  the  carbon-oxygen  core: 

y)  He4  +  C12  ♦  016  +  y  Ey  -  7.1605MeV 

NET:  4He4  016  +  E  Et  -  14.9311MeV 

The  net  energy  released  including  this  final  reaction  is 
.096772  of  the  mass  of  the  four  helium. 

If  there  is  enough  energy,  helium  will  keep  fusing  with 
heavier  and  heavier  nulcei  (016,  Ne20  ,  etc.)  until  Fe56  is 
reached.  This  has  the  lowest  binding  energy  per  nucleon  (most 
negative)  of  any  nucleus,  and  thus  is  an  ultimately  stable  con¬ 
figuration,  the  end  of  the  stellar  line. 
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APPENDIX  II 
COMPUTER  PROGRAMS 

This  appendix  presents  the  substantial  computer  programs 
written  by  the  author  in  the  course  of  this  project.  The  programs 
are  listed  in  the  order  In  which  they  would  be  used  In  solving  an 
unknown  system. 

1.  EPHi-.;  Generates  a  schedule  of  observable  events 

2.  YLEYN:  Takes  standard  star  observational  data  and  calculates 

Zya  and  Sg  to  be  used  In  Harris'  reduction  method. 

3.  TENSOP:  Operates  on  Zya  and  S0  to  give  Gu ,  the  reduction  co¬ 


efficient  tensor 
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EF-ET  27  Apr  83  23:  IS 

100  '“n  is  proaram  provides  a  schedule  o-  epienerices  in  eastern  stsnca 
110  'time  -or  a  binary  s/sismi  siven  some  initial  ep-iemeris  anc  psnoc 
120 

130  aRIN“  "INFIJT  S“AR ,  DA7Z-T:*S  GRCLP  (rG»Dr:-iR»f*,S)  r  AND  =  ERMD" 

1&0  INPUT  SS,  T*,  D,  H,  I",  S,  1= 

150  JaD+H/Z4-*-*/  1440+S/SS400 
130  3RINT 
170  PRINT 

ISO  PRINT  “A3 ( 22 )  J  "ES7-  EPHETERIDES  FOR  “  ;S$ 

ISO  aRI\“ 

200  3RIN“  “A3 (5)1  "GLT  OF  ECLIPSE":  ~A3(43):  ”'rGTAi_I“Y" 

210  aRIN“ 

220 

240  N*1 

250  FOR  J»J  "0  200  STEF  a/4 
2S0  A= I NT ( 2  ) 

270  3= J-A 
300  C  =  0 

310  READ  "1*,D1 

313  DATA  DEC,  31,  JAN,  31,  FE3 ,  28 
3 1 S  DA“A  TAR,  31,  APR,  30,  TAY ,  31,  3GRT ,  100 
320  C=C+D 1 

330  IF  C*>A  THEN  343 
340  GOTO  210 

343  IF  V,1$="3QRT"  “HEN  S7S 
245 

2^a  N=N+I 
250  D  =  A-C-rD  1 
280  H= INT  < 3*24 ) 

270  G1 =B*24-H 
280  Ta:N“(Ql430) 

390  52  =  G1*S0'-T 
400  S«  :.N”( 02*80 > 

410  IF  Ha >20  THEN  480 
440  :r  -03  “HEN  480 
450  RESE“ 

4G0  GC“G  590 
470 

480  IF  *CD (N, 4)<>0  “HEN  520 
300  RESE“ 

510  GCTO  580 

520  IF  M0D(N,2)*0  THEN  580 
525  GCSU8  BOO 

530  PRINT  TABC);  D  ’,  Ml  *  J  H  5  M  ’,  S  ’,  5  X  *,  Y  ’,  2 

540  RE3E“ 

550  GCTC  580 
560  GCSL3  BOO 

5S5  3rint  ”A8(38) ;d;mi*;h;m;s;,,-,,;x;y;z 
580  RESE“ 

580  NEX-  J 


AD- A 1 34  182  DETERMINATION  OF  INDIVIDUAL  TEMPERATURES  AND  */ 

LUMINOSITIES  IN  ECLIPSING  BINARY  STAR  SYSTEMS(U)  NAVAL  * 
ACADEMY  ANNAPOLIS  MD  R  M  CAMPBELL  20  JUN  83 
UNCLASSIFIED  USNA- TSPR- 122  F/G  3/2  NL 
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YLSYN  27  Apr  83  23: 17 

1  ’his  proaram*  in  conjunction  with  TENSOP*  recuces  photometric 
cata  by  the  multilinear  approach  o f  Harris  et.  al.  (PAS3*  >.f.  £2* 

:  pp.  287-518)  as  explained  in  the  text.  T'iis  proaram  specifically 
!  ta*es  photometric  data  and  constructs  the  Z  anc  S  tensors  for 
!  manipulation  by  TENSOP  in  douo  le-prec  i  s  i  on  3AS1CS.  “his  prasrasi 
!  utilises  3ASIC7  so  that  it  may  handle  third  ran*  tensors. 

U a r : a o 1 e  assianment  is  as  follows: 

Z  and  S  are  reduction  co-efficient  tensors. 

I  Y  is  a  dummy  composite  co-efficient  tensor*  expressed  as  a  scalar 
durina  each  loop  for  storaae  conservation. 

!  X  is  the  air  mass  tensor. 

•  N  is  the  weiahtina  factor*  w=XA-2.5*  never  explicitly  expressed. 

!  Cl.9  is  the  published  colour  index  tensor. 

•  C1G  is  the  observed  colour  index  tensor. 

X  is  tne  niaht  index*  from  1  to  M. 

!  J  is  t.ie  star  index*  from  1  to  N. 

1  and  ?  are  indices  of  Z  and  S  tensors. 

H  replaces  ooth  I  arfd  P  in  the  subproaram  TERRI. 

I  riles  S7AR3V  and  S“ARL'3  store  (3-V)  and  (U-3)  cata*  respectively. 

C3£N  #9v:  "STAR3V" 

OPEN  #L'3:  "STARUB" 

OPEN  #Ti:  "’ENSDAT" 

DIM  2(83*33)*  S(83)*  X ( 25  * *0 ) , Cl? ( 25 ) , Cl C ( 25  * *0 > 

a9INT  "INPUT  #  OF  NIGHTS*  #  CF  STARS" 5 
INPU-  M  *  N 
LS’  L*2*M+3 

REDIM  Z(L»L) *S(L) *X(N,M) » C I ? ( N ) ,CIO(N,M) 

3RIN“  "INPU’  1  FOR  (3-U)*  2  FOR  (U-3)"J 
INP'J’  R 
SELECT 
CASE  R* 1 

RESET  #8V 
DO  UNTIL  END  #BU 

INPUT  #BV:  J*,K**J,K,X(  J*K)  *CI9(J>  *CIQ(w*;<) 

LOOP 

CASS  R»2 

RESET  #U3 
DO  UN’IL  END  #U8 

input  #ub:  j**x«* j,k*x( j*;o *cip( j> »cio( j,k> 

-OOP 

SELECTEND 
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YLEYN  Page  2. 


1  All  cata  tensors  are  now  complete.  Proceed  in  calculation 
!  of  Z  and  S  tensors. 

FOR  1*1  TO  L 

for  ?*:  *-0  l 

LIT  Z(I»P)=0 
LET  S ( P ) *0 
-OR  TO  N 

FOR  K*1  TO  * 

CALL  TERRI  ( I » K,  J .  K  .  X(  ,  )  .CIPO  ,  Yi  ) 

CALL  TERRI  (P»MfJ»K»X(» )»CI?( ) »Y?) 

SELECT 

CASE  X(JtK)<>0 

LET  Z< :,P)*Z<I,P)+Yi*Yp/(X( J,K) A2.5> 

IF  I »L  THEN  LET  S  ( ? )  *S  <  P )  +C 1 0  (  J  » K  )  * Y p /  ( X  (  J  r  K  /  A 2 . 5 
CASE  X  (  J  »  X  )  =  0 
SELECTZND 
NEXT  X 
NEXT  J 
NEXT  a 
NEXT  I 

!  Z  and  S  tensors  are  now  complete.  Proceed  in  transferina 
!  t.iese  to  file  TENSDAT. 

SCRATCH  #TI 
FOR  1*1  TO  L 

FOR  P»1  TO  L 

?RINT#Ti:  2(1, P) 

NEXT  P 
NEX”  I 

FOR  P*1  TO  L 

PRINT#"-:  S(?) 

NEXT  ? 

CLOSE  #'I 


YLEYN  Page  3 


!  Subprogram  TE3.9I  calculates  reduction  co-efficient  tensors  Y. 

SL'3  TE3RI  <H,«,  J,K,X(  ,  )  »C:?( ) »  Y) 

CASS  H  =  L 

LET  Y*CIP( J ) A2 
CASE  H»L-1 

LET  YaX(J»K) *C I ? ( J ) 

CASE  H*L-2 

LET  Y=C1?( J ) 

CASE  H>K 

IF  H=f*+K  THEN  LET  Y»X<J»K)  ELSE  LET  Y»0 
DEFAULT 

IF  H*K  THEN  L£~  Y»1  ELSE  LET  Y«0 
SELECTEND 
SL3END 
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-2\SC?  27  Apr  S3  23: 19 

10  '  ""his  proaram  taxes  tne  Z  anc  S  tensors  createc  in  YLEY\, 

20  '  anc  calculates  t.ie  co-efficient  tensor  G*Z^-1;*Sf  wiici 
30  '  is  usee  directly  to  recuce  ooservations  of  unxnown  stars. 

40  '  The  inversion  and  accuracy  cnecxina  suoproarams  were  written 
50  '  by  Prof.  H.  Xaplan  of  the  Naval  Academy  mathematics  department 
SO  '  in  1976. 

30 

SO  '  Variable  assianment  is  as  follows! 

100  '  Z  and  S  are  the  reduction  co-efficient  tensors  founc  in  YLEYV. 
110  '  V  is  the  inverse  of  Z 

120  '  T  is  a  dummy  worKspace  in  subproaram  LATINO,  arc  later  Z*V. 

130  '  G  is  the  co-efficient  tensor  V*S. 

140  '  M  is  the  numoer  of  niahts  ooserved. 

150  '  E  is  the  absolute  deviation  of  T  from  tne  identity  tensor. 

1 30 

170  LIBRARY  "3ASICLIB***:  MAT I NO “ 

180  PILE  #i:  "TZNSDAT" 

ISO  DIM  Z ( 83 . 83 ) r  T ( 83 . 83 ) r  V ( 83 » 83 ) »  S ( 93  )  r  G ( 83  ) 

200  PRINT  "INPUT  #  GF  NIGHTS"; 

210  INPUT  M 
220  L«2*M+3 
220 

239  RESET  #1 

240  PGR  1*1  ”0  L 
250  FOR  ?*1  TO  L 
260  INPUT  #i:  Z(IfP) 

270  NEXT  ? 

290  NEXT  I 

290  FOR  PpI  TO  L 

300  INPUT  #i:  S<P> 

210  NEXT  ? 

320 

330  PRINT  “DGEST  THOU  WISH  TO  SEE  Z  AND  S"; 

340  INPUT  g* 

350  IF  8*0" VERILY"  THEN  410 
360  CALL  ” TENS PR I " !  Z(r),L 
370  FOR  P*1  TO  L 
3S0  PRINT  S(P); 

390  NEXT  P 
400  PRINT 
410 

420  CALL  "MA“INV":  Z< » ) »T( , ) ,V( , > fL 
423  PRINT 

430  PRINT  "DQEST  THGU  WISH  "0  SEE  Z(-1)B; 

440  INPUT  Q* 

4 30  IF  G*<  > "VERILY"  THEN  470 
460  CALL  "TSNSPRI":  V(f),L 
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*70 

475  °RINT 

400  PRINT  "DQEST  THCU  WISH  THE  ACCURACY  OF  Z(-l)  CONFIRMED" ' 

4£0  INPUT  at 

500  IF  Q»<>“ VERILY”  THEN  SOO 

510  CAL-  "MATGOGDn I  Z< » > »T( , ) ,V< , ) ,L,E 

520  PRINT  E;"IS  THE  SUM  OF  TOTAL  ABSOLUTE  DEVIATION  C”  Z*Z(-1>" 

530  ORINT  "  FROM  THE  IDENTITY  TENSOR." 

540  PRINT 

550  PRINT  “DOEST  THOU  WISH  TO  SEE  Z*Z(-1>" m, 

560  INPUT  Q* 

570  IF  Q*<>" VERILY"  THEN  600 
580  MAT  T»Z*V 

590  CAL-  "TENS PR I " I  T(,),L  i 

5S5  ] 

300  3RINT  i 

605  FOR  1*1  ”0  L  i 

310  G ( I ) *0 

620  FCR  P»1  H3  L  i 

630  G  < I ) =G ( I ) +V ( I r  P ) *S ( P )  I 

640  NEXT  P  | 

345  NEXT  I 

350  PRINT  "***  FIAT  G  ***"  ! 

355  PRINT  | 

360  FOR  I»1  ”0  L 
370  PRINT  G ( I ) 

680  NEXT  I 
SSO  END 
3S5 

700  SUB  ,,T,ENSPRI":  V(»)rw 

710  PRINT 

720  FOR  I»1  TO  L 

730  FOR  ?»1  TO  L 

740  PRINT  V(I,P) J 

750  NEXT  P 

760  PRINT 

770  NEXT  I  * 

780  PRINT  | 

790  SUBEND  j 

Rtady  ' 
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APPENDIX  III 


LIST  OF  PROGRAM  BINARIES 


SYSTEM 

a 

<5 

PER 

Nsax 

m 

“tot 

^TOTALITIES 

RT  And 

23-10-25 

52-56-02 

0.629 

4-39 

9.5 

9.7 

0 

TW  And 

00-02-25 

32-45-02 

4.12 

4-01 

9.0 

11.0 

0 

AA  And 

23-04-37 

47-35-40 

0.935 

4-30 

10.2 

10.5 

0 

V346  Aqu  20-09- li 

10-17-55 

1.11 

2-55 

9.0 

9.5 

1 

TZ  Boo 

15-07-32 

40-02-01 

0.297 

4-16 

10.5 

11.0 

40 

AC  Boo 

14-55-53 

46-25-49 

0.352 

4-28 

10.5 

10.9 

41 

SS  Can 

07-14-21 

73-21-51 

4.82 

4-56 

10.0 

10.6 

4 

SV  Can 

06-37-40 

82-17-14 

0.593 

4-28 

9.0 

9.1 

21 

SZ  Cam 

04-06-19 

62-17-17 

2.69 

4-51 

7.5 

7.7 

2 

TU  Cam 

05-53-26 

59-53-09 

2.93 

4-48 

5.5 

5.6 

6 

S  Cnc 

08-42-58 

19-05-43 

9.48 

3-25 

8.0 

10.1 

0 

RS  CVn 

13-09-50 

36-01-29 

4.80 

4-08 

8.0 

9.2 

3 

TV  Cas 

02-44-29 

65-39-21 

1.43 

4-54 

8.5 

8.6 

1 

TX  Cas 

02-50-54 

62-42-50 

2.93 

4-51 

9.0 

9.3 

1 

YZ  Cas 

00-44-30 

74-53-46 

4.47 

4-55 

5.5 

5.6 

0 

AR  Cas 

23-29-15 

58-27-21 

6.07 

4-46 

5.0 

5.1 

0 

DO  Cas 

02-4CH05 

60-28-53 

0.685 

4-49 

8.5 

8.6 

3 

U  Cep 

01-00-44 

81-47-04 

2.49 

4-32 

7.0 

9.8 

0 

ZZ  Cep 

22-44-28 

68-02-35 

1.72 

4-45 

9.5 

10.0 

0 

a  CrB 

15-33-58 

26-46-20 

17.36 

3-46 

2.2 

2.3 

1 

BR  Cyg 

19-40-24 

46-44-38 

1.33 

4-29 

9.5 

10.5 

0 

DK  Cyg 

21-34-20 

34-31-11 

0.471 

4-05 

10.5 

11.1 

0 

KR  Cyg 

20-08-25 

30-30-04 

0.845 

3-56 

9.0 

9.1 

2 

MR  Cyg 

21-58-18 

45-54-07 

1.68 

4-31 

8.5 

8.7 

0 

V380  Cyg 

19-50-02 

40-33-21 

12.43 

4-17 

5.5 

5.6 

0 

V382  Cyg 

20-18-09 

36-17-16 

1.89 

4-08 

9.0 

9.8 

0 

U  Del 

20-36-54 

18-13-27 

4.81 

3.22 

9.5 

12.5 

1 

DM  Del 

20-38-48 

14-22-02 

0.845 

3-10 

8.5 

8.7 

1 

SX  Dra 

18-04-19 

58-23-45 

5-17 

4-46 

10.0 

11.8 

1 

TV  Dra 

15-33-36 

63-57-49 

2.81 

4-52 

7.5 

9.8 

5 

RW  Gem 

06-00-26 

23-08-27 

2.87 

3-37 

10.0 

11.8 

2 

RY  Gem 

07-26-26 

15-41-45 

9.30 

3-14 

8.5 

10.5 

1 

TT  Her 

16-53-57 

16-51-52 

0.912 

3-18 

9.5 

9.8 

6 

AK  Her 

17-13-12 

16-22-15 

0.422 

3-17 

8.5 

8.9 

15 

VZ  Hya 

08-30-52 

-06-15-40 

2.90 

1-08 

9.0 

9.5 

1 

AM  Leo 

11-01-18 

09-59-21  . 

0.366 

2-54 

8.5 

9.1 

32 

SV  Lyn 

08-06-53 

41-51-03 

0.644 

4-20 

9.5 

9.8 

14 

TU  Mon 

07-52-29 

-02-59-51 

5.05 

1.44 

8.5 

10.9 

2 

UX  Mon 

07-58-27 

-07-27-29 

5.90 

0-58 

8.0 

8.9 

0 

RV  Oph 

17-33-46 

07-15-30 

3.69 

2-42 

9.5 

11.5 

1 

V566  Oph 

17-56-02 

04-59-16 

0.409 

2-32 

7.5 

7.9 

8 

W  Orl 

05-32-40 

-01-10-09 

1.49 

1-57 

5.0 

5.2 

1 

EE  Peg 

21-39-12 

09-06-25 

2.63 

2-50 

7.0 

7.1 

0 

RY  Per 

02-44-33 

48-04-24 

6.86 

4-32 

8.5 

10.3 

1 

ST  Per 

02-59-00 

39-07-25 

2.65 

4-14 

9.5 

12.0 

0 

100 


SYSTEM 

a 

PER 

ihax 

m 

“tat 

#  TOTALITIES 

U  Sge 

19-1S-04 

19-3?-45 

3758 

3-27 

6.5 

10.1 

1 

RW  Tau 

04-02-51 

28-04-50 

2.77 

3-50 

8.0 

12.5 

2 

S2  Tau 

04-35-39 

18-43-15 

0.416 

>24 

10.5 

11.1 

2 

X  Trt 

01-59-35 

27-48-23 

0.972 

3-49 

9.0 

11.4 

0 

W  UMa 

09-42-34 

56-01-49 

0.334 

4-43 

8.5 

9.2 

36 

AG  Vir 

12-00-11 

13-06-10 

0.643 

>05 

9.0 

9.4 

20 

LIST  OF  STANDARD  STARS 

SYSTEM 

a 

6 

m 

(B-V) 

(U-B) 

SPECTRAL  CLASS 

0  And 

00-16.3 

58-35 

4.60 

+.06 

+.02 

A2V 

U  Gas 

01-06.6 

54-53 

5.12 

.69 

+.09 

G5V 

8  Tri 

02-08.5 

34-55 

3.00 

.13 

+.08 

A5III 

U  Get 

02-44.0 

10-03 

4.25 

.31 

+.05 

F0IV 

6  Per 

03-41.7 

47-44 

3.03 

-.14 

-.52 

B5III 

134  Tau 

05-48.6 

12-39 

4.91 

-.07 

-.20 

B9IV 

S  Gen 

07-44.3 

28-04 

1.15 

1.00 

+.84 

K0III 

n  Hya 

08-42.5 

03-28 

4.31 

-.196 

-.74 

B3V 

X  UMa 

08-58.1 

48-07 

3.12 

.18 

+.06 

A7V 

40  Leo 

10-18.8 

19-34 

4.83 

.44 

-.01 

F6V 

S  Vlr 

11-49.7 

01-52 

3.63 

.54 

+.10 

F8V 

78  UMa 

13-00.4 

56-27 

4.93 

.57 

0.00 

F2V 

8  Con 

1>U.1 

27-57 

4.30 

.56 

+.05 

GOV 

n  Boo 

1>53.8 

18-29 

2.70 

.59 

+.20 

GO  IV 

a  Ora 

14-04.0 

64-27 

3.64 

-.05 

-.09 

AOIII 

a  Boo 

14-15.0 

19-19 

-0.06 

1.23 

1.26 

F2III 

e  CrB 

15-56.8 

26-56 

4.15 

1.227 

1.28 

K3III 

t  Her 

16-19.2 

46-21 

3.89 

-.155 

-.56 

B5IV 

6  Her 

17-14.3 

24-51 

3.14 

.08 

+.08 

A3  IV 

72  Her 

17-20.0 

32-31 

5.39 

.62 

+.07 

GOV 

Y  Dra 

17-56.2 

51-29 

2.22 

1.52 

1.86 

K5III 

a  Del 

20-38.7 

15-52 

3.77 

-.06 

-.23 

B9V 

1  Peg  A 

21-21.3 

19-44 

4.09 

1.10 

1.05 

Kill 

55  Peg 

23-06.2 

09-19 

4.50 

1.56 

1.81 

M2III 

Y  Cep 

23-38.6 

77-32 

3.22 

1.05 

+.92 

K1IV 
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APPENDIX  IV 

STANDARD  STAR  OBSERVATIONS 

All  of  Che  observations  of  standard  stars  made  In  the  course 
of  the  project  are  given  In  this  appendix.  The  first  page  is  (b-v) 
data,  and  the  second  is  (u-b) .  The  order  of  data  is: 

Star,  Date,  Star  Index,  Night  Index,  Air  Mass,  Published  Color  Index, 
Observed  Color  Index 


This  data  was  accessed  by  YLEYN  in  order  to  calculate  2^a  and  S0. 


STARS 

U  27 

A^r  83 

bGEM, 

31  JAN, 

1,1,  1 

nriYA , 

31  JAN, 

2,1,  1 

b  GEM  , 

8FE3, 

1,2,  1. 

dPER  , 

8FE3 , 

3,2,  1 . 

bGEM, 

17FE3 , 

1,3,  1 

dPER, 

i7=fc3. 

3,3,  1 

78UMa 

,  18FE3 

,  4,4, 

bGEM  , 

21FE3 , 

1,5,  1 

iUMa  » 

21FS3 , 

5,5,  1 

134TAU  ,  2MAR 

,  8,8  1 

b  GEM  , 

2MAR , 

1,8,  1. 

iUMA* 

2MAR , 

5,6,  1 . 

b  VIR » 

2MAR, 

7,8,  1. 

bGEM, 

15MAR, 

1,7,  1 

iUMa, 

15MAR , 

5,7,  1 

78UMa 

,  15MAR 

,  4,7, 

40LE0 

,  15MAR 

,  8,7, 

bCOM , 

15MAR, 

9,7,  1 

iUMa, 

23MAR , 

5,8,  1 

nBOO , 

23 MAR, 

10,8, 

aDRA, 

23MAR , 

11,8, 

eCrB, 

23MAR , 
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